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respect to expansion e-folds and density perturbations in the full adiabatic 
regime, thus insuring a valid solution regime. For particle physics, the model 
is nonstandard since it requires a large number of fields, > 10^. A particle 
physics/string interpretation of the model and solutions is discussed, which 
can accommodate the large field number requirement. 

PACS number(s): 98.80 Cq, 05.70.Ln, ll.lO.Wx 
^iep-ph/990il09l 
In press Nuclear Physics B 2000 



*E-mail: ab@ph.ed.ac.uk 

''^Present address, Department of Physics and Astronomy, University of Edinburgh, Edinburgh 
EH9 3JZ, United Kingdom 



1 



I. INTRODUCTION 



Inflation is a compelling resolution to the cosmological puzzles, because it explains how 
a large class of initial states evolve into a unique final state that is consistent with our ob- 
served Universe. The background cosmology that defines infiation is a Friedmann- Robertson- 
Walker (FRW) cosmology in which the scale factor R{t) has positive acceleration R{t) > 0. 
The importance of such a background cosmology for one of the cosmological puzzles, the 
horizon problem, has been appreciated for the longest time (for more details please 
see 0). Later, Guth in his pivotal paper [|] pointed out the importance of this back- 
ground cosmology, which he named infiation, to solving another cosmological puzzle, the 
flatness problem. Moreover, Guth's paper drew attention to the relevance of cosmological 
phase transitions in attaining a dynamical particle physics explanation for inflation. These 
ideas were foreshadowed by some earlier works 0. Furthermore, most of these ideas about 
cosmological phase transitions as well as suggestions about their importance in explaining 
the cosmological puzzles were initially developed by Kirzhnits and Linde The thermal 
fleld theory and its application to cosmological phase transitions was considered further by 
Weinberg ||^ and Dolan and Jackiw [0. 

These early ideas all converged into a single benchmark concept about inflationary dy- 
namics, the slow-roll scalar fleld scenario This picture postulates the existence of a 
scalar fleld 0(x, t), named the inflaton, which governs inflationary dynamics. At some early 
time, it is assumed that the energy density p and pressure density p were dominated by the 
homogeneous component of the inflaton foit), where, in addition, ipo(t) is assumed classical, 
with 

P = V{^o) + ^00 + Pr (1) 

and 

p = -V{^o) + + ^pr. (2) 

V{ipo) and <pl/2 are respectively the potential and kinetic energy of the inflaton and pr is 
a component of radiation energy density. The key observation for inflationary dynamics is 
when the potential energy dominates, ^^(v^o) > 0o/2;Pr, the equation of state from Eq. (|l]) 
and (0) becomes that of the vacuum p ^ —p, and this is the required condition in FRW 
cosmology for inflationary expansion. In the slow-roll inflationary scenario, considerations 
from observation generally require the stronger condition ^^(v^o) ^ 0o/2- This condition 
along with the assumption that the potential is slowly varying implies the equation of motion 
for ipo is flrst order in time. 

The conditions on pr lead to two different types of thermodynamic regimes of inflation. 
For Pr ^ expansion is isentropic during inflation, so that the Universe rapidly becomes 
supercooled. In this supercooled inflation regime, the termination of the inflationary period 
into a radiation dominated period occurs through a short intermediate reheating period, in 
which all the vacuum energy of the inflaton is converted into radiation energy. Alterna- 
tively, for V{ipo) = py > Pr > 0, expansion is non-isentropic during inflation, so that the 
temperature of the Universe may still be sizable. In this warm inflation regime, conversion 
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of vacuum energy into radiation energy occurs throughout the inflationary period and the 
inflationary regime smoothly terminates into a radiation dominated regime without an in- 
termediate reheating period. Warm inflation cosmology, i.e., inflationary dynamics without 
reheating, was formulated in [p!0| , [ll| . An earlier paper by Fang and the author [jl2| developed 
some underlying ideas, although the statement of that paper is general to all of inflationary 
cosmology. The demonstration that non-isentropic expansion, the background cosmology of 
warm inflation, can be realized in FRW cosmology from a plausible fleld theoretic dynamics 
was given in |]13HT 



Although the scale factor dynamics of the background cosmology in both supercooled and 
warm inflation regimes is similar, the microscopic scalar fleld dynamics in the two regimes is 
very different. In warm inflation dynamics, during the inflation period the inflaton interacts 
considerably with other flelds. These interactions permit energy exchange and this is how 
the vacuum energy of the inflaton is converted into radiation energy. As the flelds acquire 
the vacuum energy liberated by the inflaton and become excited, their reaction back on the 
inflaton damps its motion. To realize a viable inflationary regime, the inflaton must support 
the vacuum energy sufficiently long to solve the horizon/flatness problems, A^^e ~ 60. In 
warm inflation this implies the reaction of all the flelds on the inflaton must be sufficiently 
strong to overdamp its motion. Thus slow-roll motion in warm inflation is synonymous 
with overdamped motion. Furthermore, the effective dynamics of the inflaton is analogous 
to time dependent Ginzburg-Landau scalar order parameter kinetics. Such a kinetics has 



been derived from a near-thermal-equilibrium quantum field theory formulation |16|. In this 



derivation, the classical inflaton is deflned as the thermal expectation value of the inflaton 
fleld operator 0(x, t), v^(x, t) =< 0(x, t) >t- 

In supercooled inflation dynamics, the inflaton typically is modeled as noninteractive with 
other fields during the inflation period. Supercooling implies the Universe is in the ground 
state |0 > during inflation. It is assumed that the vacuum is translationally invariant. The 
classical inflaton that emerges in the slow-roll equation and energy-pressure densities, Eqs. 
(0) and (IP, is deflned as the expectation value of the inflaton fleld operator with respect 
to |0 >, (p{'x,t) =< O|0(x, t)|0 >. Several authors 0-|l9| have argued on the basis of 
saturating the momentum-position uncertainty principle that the quantum equations for </> 
go over into classical equations for a given inflaton mode of comoving momentum kc, (f)k^{t), 
once the physical momentum associated with the mode kp = kc/-R(t) crosses the Hubble 
radius 2n/kp > 1/H. 

The slow-roll scalar fleld scenario has a second aspect to it. Small fluctuations of the 
inflaton about its homogeneous component provide the initial seeds of density perturbations 
pil| , PT| ]. These density perturbations produced during inflation evolve into the classical in- 
homogeneities observed in the Universe. For warm inflation dynamics, the fluctuations of 
the inflaton are thermally induced. As such, these initial seeds of density perturbations 
already are classical upon inception. In supercooled inflation dynamics, the fluctuations of 
the inflaton arise from zero-point quantum fluctuations and so are purely quantum mechan- 
ical. In this case, it must be explained how these initial quantum fluctuations evolve into 
the classical inhomogeneities observed in the Universe. This problem often is referred to 
as the quantum-to-classical transition problem of supercooled inflation. Early resolutions 
to this problem followed the same reasoning as for the homogeneous component. Based on 
uncertainty principle arguments for the fleld amplitude and its conjugate momenta, these 
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arguments concluded that the fluctuations can be treated classically once the physical wave- 
length of a give mode is larger than the Hubble radius 



17-19 . Sasaki noted 



that this 

criteria for classical behavior is not invariant under canonical transformations. Thus unless 
a definite physical significance could be ascribed to the field amplitude and its conjugate mo- 
menta, such arguments are insufficient for explaining the classical realization of the quantum 
process. 

The complete resolution to the quantum-to-classical problem in supercooled inflationary 
dynamics has been understood to occur via a process which eliminates quantum interference 
between macroscopically distinguishable events. Such a process generally is termed decoher- 
ence [^. A common way to introduce decoherence is through an external environment with 
which the inflaton interacts. Various studies have examined decoherence along these lines 
in which external fields couple to the inflaton and even in which the short wavelength 



modes of the inflaton act to decohere the long wavelength modes 



For warm inflation, since dissipation is a fundamental aspect of the dynamics, decoher- 
ence is an automatic consequence. So warm inflation contains a good example of decoherence 
within an inflationary dynamics. For supercooled inflation, decoherence is not a natural re- 
quirement of the dynamics, but rather it must be imposed an an additional condition. Thus 
to the extent of relevance of decoherence, warm inflation realizes the full consequences of it 
both to yield a classical dynamics and to exploit the fluctuation-dissipation effects associated 
with it. 

In the basic warm inflation picture, the only requirement is that radiation energy is pro- 
duced during inflation and that the mechanism of production is via dissipative effects on the 
inflaton. Up to now, the only quantum field theory realization of this picture is when the 
radiation energy is near thermal equilibrium [16,26|. This is an interesting regime for devel- 



oping a warm inflation dynamics, since it is the best understood regime of nonequilibrium 
quantum field theory. For the most part, there is limited understanding of nonequilibrium 
quantum field theory. However near thermal equilibrium, the advantage is that the state of 
the system can be studied as a perturbation about the thermal equilibrium state. Further- 
more, if all macroscopic motion is slow relative to the relevant microscopic time scales, the 
macroscopic dynamics can be treated adiabatically. Finally, if the interactions are weak, 
perturbation theory is applicable. In such a thermalized, adiabatic, perturbative regime, a 
well formulated quantum field theory dissipative dynamics can be formulated. The founda- 
tions for this were developed by Kubo and Zubarev p8|, who basically examined the 
full dynamical consequences of the near-thermal-equilibrium fluctuation-dissipation theo- 



rem 



Quantum mechanical models that implement the fluctuation-dissipation theorem 
within a dissipative dynamics have been known for a long time , although in recent times 

They have been 



such models generically have been termed Caldeira-Leggett models [pl| ,f 
studied in relativistic models by various authors [0^ . 

Within a realistic scalar quantum field theory model, Hosoya and Sagagami |Q ini- 
tially formulated dissipative dynamics. In their formulation, the near-thermal-equilibrium 
dynamics is expressed through an expansion involving equilibrium correlation functions. Al- 
though their formulation is physically transparent, formally it is cumbersome. Subsequently 
Morikawa BH] formulated the same problem in terms of an elegant real time finite tempera- 



ture field theory formalism. These works were developed further by Lawrie and Gleiser 
and Ramos For warm inflation, the overdamped regime of the inflaton is the one of 
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interest. A realization of overdamped motion within a quantum field theory model based 
on this formulation was obtained in [|T^ . 

An application of a near-thermal-equilibrium fluctuation-dissipation dynamics for warm 
inflation first was examined from a Caldeira-Leggett type model . For a realistic quantum 
field theory model, the results of [|T^ were applied to warm inflation in |^ and a solution 
to the horizon/flatness problem was presented. 

The warm inflation solution in was based on a specific quantum field theory model 
that generally has been termed distributed mass models (DM-models) |[T6|j38[] . In this model, 
the inflaton interacts with several other fields through shifted couplings g'^i^Po ~ ^lYx^ 
g{'^Q — Mi)%l)i%l)i to bosons and fermions repectively. The mass sites Mj are distributed over 
some range. As the inflaton relaxes toward its minimum energy configuration, it decays into 
all fields that are light and coupled to it. In turn this generates an effective viscosity. In 
order to satisfy the e-fold requirement of a successful inflation, N^. > 60, overdamping must 
be very efficient. The purpose of distributing the masses Mj is to increase the interval for 
(fQ in which light particles emerge through the shifted couplings. 

Within this simple near-thermal-equilibrium quantum field theory formulation, the dis- 
tribution of mass sites appears necessary to sustain the inflaton's overdamped motion suf- 
ficiently long to satisfy the e-fold requirement. On the one hand, the DM- model may be 
regarded as an intermediate step towards realistic warm inflation models. On the other hand, 
the hierarchy of mass sites Mj in this model is suggestive of mass levels of a fundamental 
string. Based on this hypothesis, it was shown in [^] that DM-models can be obtained 



from effective supersymmetric theories. Further development of the string interpretation of 
DM-models is given in p9[ and [H0|. In this paper we study a wide range of warm inflation 



solutions for the DM-models, that extends the single case studied in [26|. The specifics of 



the extensions are clarified in the sections to follow. In addition, here the first estimates are 



given for thermally induced density perturbations [|12| in the DM-models. 

There is a second aspect to the present work, which perhaps is more important. It 
was mentioned above that the near-thermal-equilibrium quantum field theory formalism 
developed in [0^-|37| is unambiguously valid within a thermalized, adiabatic, perturbative 
regime. However in the basic formulation, the criteria for consistency are specified in terms 
of only a set of limiting inequalities (i.e., ^). As such, these criteria are not specific about 
the extent to which the inequalities must be satisfied. In light of this, a convincing proof that 
a given dynamics is consistent with this formalism is if the solution space of interest exists for 
an arbitrary degree of validity for the consistency inequalities. For warm inflation dynamics, 
the solution space of interest is the regime of observational consistency with respect to 
expansion e-folds, ~ 60, and density perturbations, 6p/p ~ 10^^ ^T|] (for a review of the 



basic observational facts please see p2|,|3|,|43[|). In this paper, we show that the consistency 
inequalities of the underlying formalism are satisfied to an arbitrary degree and alongside 
this, an observationally consistent warm inflation regime always exists. In particular, as will 
be seen, the most restrictive consistency conditions involve the adiabaticity requirements. 
In the solutions, it will emerge that the degree of adiabaticity can be controlled by one 
parameter, a, with adiabaticity improving as a — » 0. What will be shown is in this limit, an 
observationally consistent warm inflation regime with respect to Ne and 6p/p always exists. 
This result has a fundamental relevance, since it is an existence proof within quantum field 
theory for observationally consistent inflationary dynamics from start to finish: from an 
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initial radiation dominated or inflationary regime, to an inflationary regime and finally into 
a radiation dominated regime. 

The paper is organized as follows. In Sect, ^the DM model Lagrangian is presented and 
useful results about its effective potential are reviewed. In Sect. |ITl| the basic equations of 
warm inflation and the consistency conditions on the solutions are reviewed. Subsect. |1I1 C 
provides a convenient summary of all parameters, notation, and terminology used in this 
paper. In Sects. ^ and |V| the solutions for e-folds and density perturbations respectively are 
presented. Subsects. |IV C| and |V C| discuss general features about the solutions in order to 
ease the effort to understand the calculations. This paper will not focus on phenomenological 
consequences of the model. In Sect. [V^some example applications of the solutions are given. 
Subsect [VI A| computes the dimensional scales of the relevant quantities in the cosmology. 
Subsect. |VIlj| analyzes this warm inflation model in the limit of arbitrary adiabaticity. 
Subsect. |VI C| examines the solution's dependence on A, the self coupling parameter. In 
the concluding section |V11| , first improvements to this calculation are discussed. Second, 
a particle/string physics interpretation of the model is discussed. Finally some concluding 
perspectives are given about the model and solutions. 



II. MODEL 

Consider the following Lagrangian of a scalar field interacting with Nm x scalar 
fields Xik and Nm x fermion fields ipik, 

— — 1 77?^ A 

£[0, x^k. ^ik. ^iu. X\. = 2 - -Y<P^ - - H - VM) 

«max r 1 f „2 ^ 

+ E E \i^x^uf - fxt, - f (0 - M.f xl 

% — ^min. ^ — ^ ^ 
«max ^ii _ 

^ — ^min k = \ 

+ E [-^{d,^f - + E {#1 • (3) 



This model is described in [0,0 and named the distributed mass model (DM model), since 
the interaction between with the Xik and ipik fields establishes a mass scale distribution 
for the Xik and ipik fields, which is determined by the mass parameters {Mi]. 

The finite temperature effective potential for this model was computed in 0,^] for a 
nonzero homogeneous field amplitude (0(x, if:))^ = ^oif). Because of the shifted coupling 
arrangement, the self-coupling parameter A is not corrected by x ~ and "0 — interactions. 
However, these interactions will generate new terms oc (0 — Mj)" and they are discussed in 
the next section. The finite temperature, field dependent mass are 

ml^{^o.T)=g\^o-M,f + ^^l^^{T), (4) 
mJ^^((^o, T) = [g{^o - M,) + ^i^,XT)f (5) 
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and 



mJ(^o,T) = m^ + M + ^J(T). (6) 



/i(T) are thermal mass corrections with ^.^{T) ~ \/~XT and fi^^^{T) ~ gT. These thermal 
masses are nonzero only when the field dependent contribution to the respective particle's 
mass is below the temperature scale. 

Note that the Xik and ipik effective field-dependent masses, m^^^,^^.i^{(po,T), can be con- 
strained even when (0) = ip^ is large. The (px, (pip interactions can be made reflection 
symmetric, (f) —(p, but for our purposes we consider all Mj > and y^o > 0. The parame- 
ters Mi will be referred to as mass sites. The x and ip flelds will be referred to as dissipative 
heat bath flelds. 

For our calculations we will consider all coupling constants to be positive: A, 
fik,gfk^hik, fi > 0, and for simplicity choose them to be f- = fit = /, gik = hik = g. 
Also, we will set A^^^ = N, = N/A, which implies an equal number of bose and fermi 
degrees of freedom at each mass site Mj. This relation on the degrees of freedom along 
with our choice of coupling implies a cancellation of radiatively generated vacuum energy 
corrections in the effective potential For convenience we will choose the mass scales 



to be evenly spaced as Mj = iM/g, i = Vin, • • • , Vax = Vin + Nm, where M is the mass 
splitting scale between adjacent sites. Here Mj^.^ and Mi^^^^ bound the interval for ip^ in 
which dissipative dynamics, thus warm inflation, is possible. 

The Lagrangian Eq. @ also contains Nr bosonic x"^ and Nr/4: fermionic ip^ flelds that 
do not interact with the inflaton (p. These flelds represent additional degrees of freedom 
in the heat bath that otherwise do not contribute to the dissipative dynamics of (p. These 
flelds will be referred to as non-dissipative heat bath flelds. Our choice of equal numbers of 
bose and fermi degrees of freedom for the non-dissipative heat bath flelds is not necessary, 
but done here for notational convenience. 

The non-dissipative heat bath flelds as written in the above Lagrangian are completely 
noninteracting with either (p and the dissipative heat bath flelds coupled to (p. Thus, in 
the form written, none of the energy liberated by the inflaton (p is transferred into the 
non-dissipative heat bath flelds. To properly realize such heat bath flelds, decay channels 
are necessary between them and the system of flelds coupled to (p. We will not attempt to 
correct this shortcoming but assume that it is possible. Up to the level of approximation in 
P7| , P!B|P^ , which is what is applied here, this problem can be solved. Here, simply we want 
to explore the consequences of this possibility. As will be seen, observationally consistent 
warm inflations can be achieved with no non-dissipative heat bath flelds A,. = 0, and in many 
cases this is the optimal situation. However, as also will be seen, interesting possibilities 
arise when there are a very large number of non-dissipative heat bath flelds. 

The Lagrangian Eq. has two extensions to the one in [^]. First the above Lagrangian 



allows for an overall constant shift in the vacuum energy Vq. Second, an additional term has 
been added to the interaction potential, Vi{(p). This additional term is deflned to be zero 
within the interval of (po in which we study warm inflation, Mj^.^ < (fo < Mi^^^. However 
outside this region, Vi{(p) will be suitable to permit an absolute minimum of the <^-efffective 
potential at zero potential energy. In the calculations to follow, one case is treated where 
Vq and Vi(0) are nonexistent, Vq = and Vi{(p) = for all (p. 
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A few comments are in order to motivate the extensions Vq and Vi(0). The focus of this 
study is to understand the quantum field theory dynamics that underhes warm inflation. 
In this respect, the extensions Vq and Vi(0), permit exploration of warm inflation dynamics 
in a variety of interesting regimes. In this paper, we will not attempt to motivate these 
extensions from particle physics. However, note that the effect of Vq and Vi(0) is similar 
to a plateau region that could arise from higher polynomial or non-polynomial potentials. 
Non-polynomial potentials can occur in particle physics from nonperturbative effects. For 
example, nonperturbative mechanisms for dynamical SUSY breaking involving instantons 
can yield non-polynomial potentials (for a review please see [Q). Finally note that the 
basic results in this paper do not rely on these extensions Vq and Vi(0). A warm inflation 
regime consistent with quantum field theory and observational requirements on e-folds and 
density perturbations can be obtained without Vq and Vi(0), i.e. Vq = Vi(0) = 0. This case 
is treated in Subsects. [IVB 2| and |VB2 . 



III. STATEMENT OF THE PROBLEM 

Let t = tsi = define the time when warm inflation begins (BI) and t = Iei define 
the time when warm inflation ends {EI). The thermodynamic state of the Universe during 
warm inflation is given by the time dependent temperature T{t). The results presented in 
the next sections require two temperature scales Tbi = T(tBi = 0) and Tei = T(tEi), the 
temperatures of the universe at the beginning and end of warm inflation respectively. These 
two temperatures will be expressed through the ratios 

J-BI 

and 

KM = ^, (8) 

where M, defined in Sect. |I|, is the splitting scale between adjacent mass sites in the DM 
model. 



A. Basic Equations 

The basic equation for warm inflation dynamics is an effective equation of motion for 
the scalar field cpo, which in the warm inflation regime describes overdamped motion. This 



equation has been derived for DM-models in with details about the bosonic disspative 



term in and the fermionic dissipative term in |^ . Intuitive explanations of the formalism 
are given in |3^,^^. In our calculations, the contribution from fermionic dissipation is 
ignored. The only effect considered from fermions is their contributions to the temperature 
dependent yj-effective potential. The most important term they contribute is their T = 
radiatively generated vacuum energy term, since, due to our choice of bose-fermi degrees 
of freedom, it cancels a similar term from the bose sector. The high-T modifications from 
the fermion fields are not essential for the overall consistency of the calculation and they 
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also do not cause any new problems. We will account for these contributions in the ip- 
effective potential. Had we accounted for the fermion effects to dissipation, the damping of 
(fQ would increase, which in turn would increase the robustness of warm inflation solutions. 
However in our case, where there are one-forth the number of fermion versus bosonic fields, 
for the most part, based on the high-T expressions for the dissipative terms p!6|j2^j45[1 , their 



contribution would enhance disspation by only 20 percent. We believe one can configure 
situations in which the fermion versus bose contribution to dissipation dominates, which 
is interesting since the functional behavior of the two types of dissipation differ. This is 
another possibility to examine, especially since the analysis in found greater success 
with fermionic dissipation. However here we will not pursue that direction. 



For all the cases studied below, the basic v?o effective equation of motion is from |116| , |26 



where t.e. means sum over all thermally excited sites. Everything multiplying d(po/dt will 
be referred to as the dissipative coefficient and ri^{T) will be referred to as the dissipative 
function. Here V{(po,T) is the effective potential for ip^. Both the dissipative coefficient and 
(/^-effective potential emerge upon integrating out the dissipative heat bath fields Xik and 
ipik. For Xik - fields with m^. < T, the high-T expression for rjf^ is 

- 1'^ (10) 



7iTC{gjy 



where 

Cig,f)^g' 



rLi2 



7r2 \m^{ipo,T) 



£2 £2 

+ ^^^?' + ^ (11) 



and 7]' = 48 ln(2T/yU^). In an expanding background, the interaction of (p with the metric 
also yields a SHipo term in the (po-eSective equation of motion. As discussed in [T^J2^, this 
term is dropped in Eq. (j^) because the thermal dissipation term dominates. 

In order to be within an inflationary regime, V{ipo, T) must be vacuum energy dominated, 
so that V{ipo,T) pa Pi,((y9o), where Pvifo) is the vacuum energy of the 0-field, Pv{,^o) = Vq + 
P^ifo), where Vq is the vacuum energy density shift parameter and p^ifo) = V{fo,T = 0). 
In this case, the rate of decay of vacuum energy is obtained as dp^/dt = V'{ipo)ipo. The 
energy balance equation describing the transfer of vacuum energy to radiation energy is 

Pr = -AHpr - p^. (12) 



In the warm inflation regime pr decreases slowly, so that to a good approximation Eq. (p!2D 
becomes 

AHpr ^ -p.. (13) 



For the four cases considered in Sects. |I^ and 0, we examine the regime where the A^Jq 
term dominates the effective potential. 
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B. Consistency Conditions 



The validity of the basic equations of warm inflation given above requires four consistency 
checks, which we refer to as the thermalization, adiabatic, force and infra-red conditions. 
Below, these conditions are explained and exactly specified. 

The thermalization and adiabatic conditions both address the consistency of the macro- 
scopic warm inflation equations with the underlying microscopic quantum field theory dy- 
namics. Within the simple dissipative quantum field theory formulation in |T6|, the decay 
widths r,^, r^, and characterize the time scale of microscopic dynamics. The basic self- 
consistency requirement is that all microscopic time scales are faster than all macroscopic 
time scales. Thus the slowest microscopic time scale is the important one for checking self- 
consistency. For dissipative dynamics, F,^ does not play a significant role. Since we are 
restricting to the bosonic contribution to dissipation F^ is the relevant microscopic rate. 
From \TB\, the high temperature expression for it is 



rx(T) ^ —CigJl (14) 

where C is defined in Eq. ([Tl|) and F' = 1/192. For x-fields at site i, their decay rate is Eq. 



(14) when m^.^, ~ T and negligible when m^^^, ~ T. 
Hereafter, rescaled time is used 

r ^ F^(M)t (15) 

The Hubble parameter at the beginning of warm inflation will be expressed as 

Hbi = apr^iTBi) = aPKMT^iM), (16) 

where the additional parameter a is introduced and /5 and km are defined in Eqs. (^ and 
(|) respectively. 

Warm inflation fundamentally has two macroscopic time scales, the Hubble expansion 
rate and the rate of change of the inflaton, which expressed in terms of rescaled time Eq. 
(p!5|) is 00 = ^ xi^y^o / ^'^) ■ ^^^f consistency requires the thermalization condition 

Hir) « F^(T) (17) 

and the y^o-adiabatic condition 



Since the fermion contribution to dissipation is being ignored, the behavior of is not very 
essential. If < F^, then the tp-iields may not be thermally excited, which in turn would slightly 
modify the temperature dependent terms in the effective potential. However since the ■0-helds 
are free from the crucial self-consistency requirements for producing dissipative effects on ipo, if 
< Ty^, is is easy to increase their number of decay channels, which thereby can increase F^ to 
the level of F^. 
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1 (^g) 



ipo dT r^(M)' 

For the cases in Sects. and|V|, H{r) changes httle (less than a factor 5) from Tbi to Tei, so 
that Hbi is its approximate scale throughout warm inflation. Therefore, the thermalization 
condition Eq. ([T7| ) requires a < 1. 



A stronger adiabatic condition also requires that the (pQ dependence in all Boltzmann 
factors varies slowly relative to the thermalization rate. The relevant thermal excitations in 
this model are the Xki and ipki fields when m^|,.^^|.^ < T, for which the stringest form of the 
thermal-adiabatic condition is 

T dr T r^(M)- ^ ^ 

The estimates to follow are based on this criteria, but it is much more restrictive than 



necessary for two reasons. First thermal excitation generally initiates once m^^^ ~ lOT. For 
the DM-model this implies mass sites have quite a long time to excite before their important 
role arises, which is when m^^^ ~ T at which point they drive the viscosity. Second, within 
the thermally excited regime m^^^ ~ T, the Boltzmann factor is saturated at ~ 1, so is 
insensitive to mass variations. 

The next consistency check is the force condition. In our case this requires the 
term to dominate the effective equation of motion. Due to vacuum energy cancella- 
tions between bosons and fermions, the one-loop T = radiative corrections [fzl, ^ , ^ 



~ Ef g^Ni^o - M,)^ Ef g^Ni^o - M^f ln[((^o - M,)V/iLorm], are eliminated where 
the sum here is over all mass sites.. The T-dependent force terms from the effective po- 
tential are 1. {g^N /S)T^Y.YiVo - M,), 2. {g^N / {ATi))TY.Tiv^ - ^^)l<^o - M,|, and 3. 
[g^N / {QA'k))Yj\^'{,Vo — \og{T / n^enovm) ■ Our approximation makes a sharp division be- 
tween thermally excited and unexcited sites. If the mass of a site fn'^.^.^tp.,. ~ g'^ifo — Mi)'^ < 
T^, that site is thermally excited and for > T"^ it is thermally unexcited or cold. 

Here we assume (po is always surrounded by sufficiently many mass sites on both sides so 
that there are no edge effects. 

There are two general features about all three types of force terms. First, the direction of 
the force flips for mass sites on opposite sides of ipo and this implies considerable cancellation 
in the summations. Second, all three forces are periodic in i^q with periodicity M. As ipo 
traverses any M-interval, it experiences an identical force profile in both directions. Thus 
the average force experienced by (po from all three force terms is zero after traversing any 
M-interval. 

In our estimates, we use the most stringent consistency condition, which is to estimate the 
largest effect from these force terms and demand that the X(Po/Q term dominates it. We will 
find a large warm inflation solution space with this stringent condition. However estimates 
based on this force condition are overly conservative, since realistically <^o is not sharp, as 
we are treating it, but rather is smeared over some interval AipQ. Smearing would imply the 
average force experienced by the v^o-packet diminishes due to the directionality dependence 
of the force terms. In fact if A(fQ > M, the force experienced by the wavepacket from these 
three force terms would be negligible. Nevertheless, for simplicity of the calculation, here 
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this fact will not be treated and we will estimate the the upper bound of the three force 
terms for a sharp lpq. 

Amongst the three force terms, note that 1 is the largest. Ignoring the premultiplying 
factors (which also is largest for 1), since for all thermally excited sites (7|<y9o — Mj| < T, it 
implies Ef (<^o - Mi)T^ > 9{^o - M,)\ipo - M,\T > g^^o - Mi)\ Thus provided 
A(/3q/6 dominates force term 1, it dominates them all0. Approximating 

gY:{'Po-M,)^T (20) 

i 

we obtain the force condition 

^>^i:(.o-M.r.fr' (21) 

The final consistency check is the infra-red condition. For the Minkowski space quantum 
field theory formalism used here to be valid, the Compton wavelength of all particle exci- 
tations should be much smaller than the Hubble radius 1/H{t) during inflation. For the 
Xiky'ipik fields, this condition is easily satisfied since their masses are generally m^.^^^-^ ~ T 
and T ^ H. This follows since for most of the time the v9o(''")-induced portion of their 
masses is large. However, even when (/?o(t) ~ Mi, the thermal mass contribution ~ gT 
generally is much larger than H. On the other hand, the 0-mass is more concerning, since 
the self-coupling parameter A is usually tiny in warm inflation. Thus the essential constraint 
that must be checked to enforce the infra-red condition is 

m^{ipo,T):^H. (22) 

In total, there are five consistency conditions. The thermal-adiabatic, <y9o-adiabatic and 
thermalization conditions all are adiabatic conditions. The latter, equivalently stated, re- 
quires the cosmological expansion rate to be adiabatic relative to the particle production 
rate. The force condition is not fundamentally required, but we have imposed it to simplify 
the ipo equation of motion. 

C. Summary of the Parameters 

In Sects. 1^ and two cosmological parameters will be computed, the number of e- 
folds Ne and amplitude of density perturbations Sp/p, in terms of the microscopic and 
thermodynamic parameters of the model and one overall scale which will be the Planck 
scale Mp. Since the calculations are fairly detailed, for convenience here all definitions of 
parameters and related terminology are summarized. 

The following terminology is used in this paper: 



^When all bosons and fermions are thermally excited, due to our choice of coupling and of ratio 
between bose and fermi fields, term 3, in fact, cancels amongst the bosons and fermions [^. 
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dissipative heat bath - fields Xifc, "^Pik that are part of the heat reservoir and 
fields have significant effect on the dissipative dynamics of the 

inflaton. 

non-dissipative heat - fields Xi ? V'l t^^^ only are part of the heat reservoir and 
bath fields have no significant effect on the dissipative dynamics of 

the inflaton. 

mass site - "location" of dissipative heat bath flelds in the DM model 

specifled by the mass parameter Mj = Mi/g, with i de- 
noting the number of the site. 

dissipative coefficient - the factor multiplying dipo/dt in the (po equation of mo- 

tion, Er-(^o-M,)V^(T) 

dissipative function - rjf{T) = rj' g'^ / (tcTC) with rj' and C deflned below Eq. 

m 

The microscopic parameters have all been deflned in detail in Sect. 0. Briefly, they are 
as follows: 

N - number of Xik flelds at every mass site Mj. Correspond- 

ingly the number of ipik flelds at every mass site is N/4. 

Nr - number of nondissipative heat bath flelds 

Nm - total number of mass sites crossed by the inflaton ipo 

during warm inflation. Ignoring small corrections at the 
two end points, this is equivalent to the total number 
of mass sites that were thermally excited at some point 
during warm inflation. 

A - (j) self coupling parameter 

g - 4> - X coupling ~ g'^/2 , 4> - ip coupling ~ g 

f - X self coupling parameter 

Vq - vacuum energy density shift parameter 

M - splitting scale between adjacent mass sites g\Mi+i — 

Mil = M 

m^^^^^^^ - mass of Xik or ipik flelds given in Eqs. (|), (||) 

- mass given in Eq. (||) 

The thermodynamics of the warm inflation is expressed through the following quantities: 

Tbi - temperature at the beginning of warm inflation 

Tei - temperature at the end of warm inflation 

hm - Tbi/M 

f3 - Tei/Tbi 

The cosmology is expressed through the following quantities: 
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R{t) - scale factor, with -R(O) = 1 

Hbi - Hubble parameter at the beginning of warm inflation, 

Hbi = H{tbi)- In the four cases examined in this pa- 
per, the scale of H{t) is of order Hbi throughout the 
inflationary period. 

a - Hbi = C(l3K,M^chi{M). It will be seen that a is the adia- 

batic parameter with adiabaticity increasing as a ^ 0. 

Pv{t) - vacuum energy density = Vq + Ptpij) where is the v?o 

dependent portion. 

Pt{t) - radiation energy density 

Ne - number of e-folds 

5pl p - amphtude of scalar density perturbation 



Some additional notation used in this paper is as follows: 

BI - begin warm inflation 

EI - end warm inflation 

r^(M) - decay width for x-helds when they are thermally excited, 

explicit expression in Eq. (p^) 
T - rescaled time r = r^(M)t 

rit.e. - number of thermally excited (i.e.) mass sites 

0(x, r) - quantum inflaton field operator 

V9(x, r) - classical inflaton field (0(x, r))^ = (y9(x, r) = v^o(t) + 

(5(y9(x, r), where v^o(t) is the homogeneous background 
field 

y - slope parameter for ip^. Since the evolution is over- 

damped, in all cases dipQ/dr oc yipo. 

kp - freeze-out momentum for density perturbations, defined 

Eq. dlli 

^^'h{'^) ~ amplitude of scalar field fluctuations, defined Eqs. ( |122D 

and ([T23|) 

kp, kc - physical {p) and comoving (c) wavenumber. By our con- 

vention they are the same at the end of warm inflation 
tei, kp(r) = kc/-R(r - tei). 

^V'k.X'^) - Fourier mode of the inflaton field (5v9(x, r); equivalently 

this mode may be expressed as S(p(kp,T), where the re- 
lation between kp(r) and k^ is specified. 

Finally, for the following expressions, we simply state where they are defined: rj' = 
481n(2T/m^) ^ 48 - below Eq. (|lTD; F' ^ 1/192 below Eq. (|14D; C - Eq. (|Tl|), also see first 
part of Sect. |VI|; T - Eq.([n5|). 
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IV. SOLUTIONS FOR E-FOLDS Ne 



In this section four warm inflation solutions are presented. The solutions we examine 
are for cases where Tbi 3> Tei in Subsect. [IV A| and Tbi ~ Tei in Subsect. |1V B| . We seek 
solutions in which N^. e-folds of inflation occur with the temperature of the Universe Tbi 
and Tei respectively at the beginning and end of warm inflation. The temperatures Tbi and 
Tei are parameterized by (3 and km in Eqs. (0) and (||) respectively These conditions on 
the solutions imply two parameters of the model are determined in terms of the others. We 
choose these two quantities to be A and ipBi/M. Under these conditions the solutions for 
'■Po{t) and T{t) are given as well as expressions for A and (pBi/M. In addition, for ease of 
reference, expressions are written for m'^/T'^, p^^r)/ pr^r), dipo/dr and T^{t)/T^{M). Finally 
the consistency conditions from Subsect. [Ill B| are evaluated. The calculation was designed 
so that the thermalization condition Eq. (|1^ is satisfied simply by requiring a < 1. The 
adiabatic, force, and infra-red conditions result in two essential inequalities that restrict the 
self-consistent region of the parameter space. No overall scale is chosen in this section, since 
the results leave this choice completely arbitrary. 

Although all the calculations are very simple, due to the generality of the results, the 
final expressions may not appear transparent. It must be appreciated that despite the 
compactness of the final expressions, they contain the solutions under the very general 
situation in which the parameters can be varied over a wide range and in which four scales 
are adjustable, the initial temperature, the final temperature, the Hubble expansion rate 
and the duration of warm inflation. 



A. Tbi > Tei 

In this regime when the temperature from the beginning Tbi to the end Tei of warm 
inflation changes significantly, the time dependence of the number of thermally excited sites 
is treated in the y^o-equation of motion. One general parametric constraint for this regime 
is 

/3 < 1. (23) 

At temperature T(r), any Xik fields with mass tti^,^ ~ g'^i'^o ~ M-if' < T"^ is thermally 
excited. Due to our choice of spacings between mass sites Mj, this implies approximately 
T{t)/M sites adjacent to the field amplitude v?o on either side are thermally excited. In our 
approximation for all sites with m^^-^ > T, the contribution to the dissipative dynamics is 
ignored. Eq. (0) therefore becomes 



vr^ dr 6 

Also, we allow for an additional bosonic and N^/A fermionic non-dissipative heat bath 
fields that contribute to the radiation energy density but have no affect on the dissipative 
dynamics of the inflaton. Thus the radiation energy density at any instant during warm 
inflation is 
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Pr{r) = ^^[N^ + ^-^]T\r). (25) 

We examine the solutions for the cases A^,- = and Nr ^ 2Tbi/Tei in Subsects. |1V A 1| 
and [IV A2| respectively. For both cases, the vacuum energy density shift parameter will 



be Vo = \ip%i/24:, so that the vacuum energy is approximately constant throughout warm 
inflation. In this case the Hubble parameter, H = ^/8^TGp ~ y/8nGpv, changes by an 
0(l)-factor during warm inflation and the number of e-folds is 

Ne ^ H Bit EI = af^KMTEI- (26) 



1. Nr = 

In this case there are no non-dissipative heat bath fields. 

Solutions 

Eq. (H) and Eq. (H) imply 



Substituting into Eq. (p4D, we find the solution 

Mr) ^ j^^, (28) 

where y = [7r222/7(33A7(r/'r)^/')] [(a/3KM)'AV(Ar3(^iO)]i/^(v9B//M)2/7. From Eq. (|2B|) for 
A^e and Eq. (0) for /5, 

,= ^^^(1-/?V3). (29) 

By equating this expression for y with the one in terms of the parameters in the model, 
and applying Eq. (|^) to eq. (pTf), two parameters of the model are determined, of which we 
choose 

) A^l 

and 



^BI 



,s„r)./» Ml 4f . h'^] (31) 



Using Eqs. (|0D and (gj), Eq. <^ becomes 
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Tir) ^ flH^V" ( ^'^^^'^^^i^-m'^ ] MVv«/V). (32, 



27r' 



Eqs. ( p8D and (|32| ) are the general solutions. 

Based on this solution, next some useful expressions are given. The number of sites that 
are thermally excited at a given instance is rite = 2T{t)/M = 2KM/{yT + 1)^. The number 
of mass sites ipQ crosses during the warm inflation period is 



with 



The mass of is 



^<^- - 1^ 2712 j ;v,(yr + 1) • ^^^^ 



The ratio of vacuum to radiation energy density is 



M^= 7(1-/?V3) (^^+^) (36) 



T7^(?/^ + l) • (37) 



p,(r) 7(1 -/3V3) 
Finally, for examining the adiabatic condition 

rf<^o 7a/52/3^^(i _ ^1/3) / \ aV2^V2^^ ^(^ 



and 

rx(T) _ /tM 

r^(M) " (yr + l)3- 



(38) 



(39) 



Consistency Conditions 

The above solutions still are subject to consistency checks. For the yjo-adiabatic condition 
Eq. (|18[) , since the thermalization rate Eq. (|39| ) decreases faster than {dipQij) / dr) / ip^ from 
Eqs. ( |28| ) and (|38D, the most stringent test is at tei'-, if the y^o-adiabatic condition holds at 
tei, then it is satisfied better at earlier times. We find the <y9o-adiabatic condition Eq. (p!8D 
requires 7a/{2Ne) < 1. Since a < 1, for this to hold, it is sufficient that 
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N, > 7/2. (40) 

The thermal-adiabatic condition Eq. (|T^ implies from Eqs. (P^, (^) and that the 
most stringent time is Tei with the constraint 



< 1. (41) 



By substituting for rj' and F', it gives a < 1/199 (at r = the condition imphes an additional 
factor of inEq. (|T]) ). 

The force condition Eq. (|2T|) imphes 



4^2r]'ra^/^pKM > 1. (42) 
Finally the infra-red condition Eq. ( P^ implies 

2lV/?(l-/5V3)^^^iV(72 



2F' «iV,((?4 + j2/8)^ 



> 1. (43) 



For the force and infra-red conditions, they also are evaluated at the most stringent instant 
during warm inflation, which turns out to be again r = tei- 

2. Nr > 2NT/Tei 



In this case the Nr term dominates the radiation energy density in Eq. (p5|). 

Solutions 

The procedure for solving this case is similar to the above case. The results are 

= ipBiexp{-yT) (44) 

and 

= 3V3(^VoVa [aP.MNrNg^] ^r) = >^mM eM-yr), (45) 

where in terms of the parameters of the model y = [tt^/ {3^/^2{ri'T')'^/^)] [XNra(3KM / {N'^g'^)Y^^ 
and in terms of the expansion e-folds from Eq. (^) 

a/?KAf ln(l//3) 

y = j^^ ■ (46) 

From the specified conditions on temperature expressed through a, /?, and km, two param- 
eters in the model are determined 

^~[~^^) NfNr ^ ^ 
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and 

M \ A2 J 2(r7T')i/2ai/2/?i/2in(l//?)Ari/2(^- ^ ^ 

Based on these solutions, some useful expressions are as follows. The number of mass sites 
that are thermally excited at a given instance are nt,e. = 2^^/ exp(— yr). The number of 
mass sites that ipo crosses during the warm inflation period is 



with 



TT 



2 



™" 2(r/T')i/2ai/2ln(l//5)ArV2- 
The mass of the (p particles is 



2 Xifl 3r]'T' af3KMH^/(3)Ng' 



The ratio of vacuum to radiation energy density is 

(t) 1 



Pr{T) (r7T')V2 1n(l//3) 
Vn 1 K 



p,(r) (Vr')i/2 1n(l//?) 
For examining the adiabatic conditions 



r^(M) 



(50) 



(52) 



exp(%r). (53) 



^ = = 2(vfO^ A^V^^ ^^^^ ^^^^ 

and 



fi;Mexp(-yr). (55) 



Consistency Conditions 

The parametric constraints from the consistency conditions are the (/^o-adiabatic condi- 
tion Eq. ([TsD with sufficiency condition 



iVe > 1, (56) 
the thermal-adiabatic condition Eq. (|I9|) 
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the force condition Eq. (^Tj) 



1/2 



< 1, (57) 



Kirf'^i^^T >l. (58) 



> 1. (59) 



and the infrared condition Eq. (|2^) 

r aN,{g* + /V8)2 

It should also be recalled that the basic requirement for this regime, A^^-dominated, is 

^ > 2fi:M. (60) 



B. Tbi ~ Tei 

For this regime, we adopt the criteria 

(3 ~ 0.5. (61) 

In this regime the number of thermally excited sites in the yjQ-equation of motion is approx- 
imated as constant, so that in Eq. (D Ei'^'(¥'o - Mifrif^iT) ^ kI^M^Nt]' /{g^T). All other 
approximations are the same as in the previous section. Thus the y^Q-equation of motion is 

iT'Kl,g'NM^ dvo _ X^l 
tt'T dr ~ 6 ■ 

Both cases studied below are for arbitrary number of dissipative and non-dissipative heat 
bath fields so that the radiation energy density is 

p,,(r) = ^{Nr + 2Nkm)T\t). (63) 

The two cases examined in order are Vq > and Vq = 0. The latter case is the DM-model 
examined in [Q. However here the calculation is extended to permit an arbitrary mass 
splitting scale M (in [^] the mass splitting scale was restricted to M ~ T) and to treat the 



time dependence of the temperature in the yjo-equation of motion. 
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1. Vb > 



Similar to the previous subsection, we choose Vq ~ A(y9^j/24. The number of e-folds 
follows from the relation Eq. (p6|). 

Solutions 

The procedure for solving this case is similar to the above cases. The results are 
and 

y> {97^'rYr^\apn\,{Nr + 2KMN)Ng^) M {yr+iy/^' ^ ' 

where in terms of the parameters of the model y = [27i'^ / {3^^^{r]'T')^^^][\^ /{aPK^^^{Nr + 
2K,MN)N^g^Y^^{(fBi/M)^ and in terms of the expansion e-folds from Eq. 



(66) 

From the specified conditions of temperature expressed through a, P, and km two parameters 
in the model are determined 

\ 87r6 J (3^N^{Nr + 2KMN) ^ ^ 

and 

^ = rq ( (y'P'^''MiNr + 2KMN)Ng^ y^' ^ 2ti^ (^-^'^t^jj^N.^Nr + 2KMNfl^ 

M ' \ A2 ) (r/T')i/2 aV2(i_ ^2)^1/2^ 

(68) 

Based on these solutions, some useful expressions are as follows. The number of mass sites 
that are thermally excited at a given instance are nt.e. = 2/tM/(?/r + 1)^/^. The umber of 
mass sites that ifo crosses during the warm inflation period is 

^ 27r^ py\l-py')4'N^{Nr + 2^MNY/' 

^ (r/T')i/2 ai/2(i_^2)jvi/2 ' ^^^> 



with 



. ^ 27r^ (3^K^i^N,{Nr + 2KMNYl^ 

{ri'Vy/^ ai/2(i_;52)jvi/2 " ^'^> 



The mass of the (p particles is 
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The ratio of vacuum to radiation energy density is 

(r) 4 Z^^iVe 



Pr(r) 34/3(r^T')2/3 (1 - 



(yr+1) (72) 



{yr + iy. (73) 



p,(r) 34/3(,^/p)2/3(l_/52) 

For examining the adiabatic conditions 

dip, a{l-l3^)KM Mr) a^l^ ^^'^ K^i^ {Nr + 2KMNf' 



2 



and 

rx(T) _ 

r^(M) " (yr+ 1)1/2- 



■T(r) (74) 



(75) 



Consistency Conditions 

The parametric constraints from the consistency conditions are the v^g-adiabatic condi- 
tion Eq. ([TSD with sufficiency condition 

Ne>l (76) 



the thermal-adiabatic condition Eq. (|T9|) 



t/?(iV, + 2/€MiV) y^' ^ ^ 



(77) 



4(^T')^/' |^ ^/^^M(iV^+2^MiV) ^ ^ ^^g^ 



2(vr')i/2 \^ /.^^jv ; 

the force condition Eq. (|2T1) 



and the infrared condition Eq. ([2^) 

377' (l-/3')/tMiV^7' 



4r'a/32iV,((74 + /2/j 



> 1. (79) 
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2. Vb = 



In this case, there is no shift parameter Vq nor extra function Vi((/)) in the Lagrangian 
Eq. (^. This is the model examined in [^], except here the calculation is extended to treat 



the time dependence of the temperature in the dissipative function. Also only examined 
the regime where the mass level splitting M ~ T, whereas here the relation between M and 
T can be varied through the parameter km- 

There are two differences in this calculation's procedures compared to the previous three 
cases. Both differences arise because here the Hubble parameter is ^^{t) dependent 

E{t) = a(3KMT^iM)^^4^. (80) 

V^BI 

First this dependence must be treated in the energy conservation equation ([T3|). Second the 
scale factor no longer grows exactly exponentially. However this calculation is also based on 
the assumption that the temperature during warm inflation does not change significantly, 
P ~ 0.5, and as will be seen this also implies ipoir), thus Pv{t), does not change significantly. 
In this case, the scale factor grows quasi-exponentially with e-folds 

r^Ei fTEi ojnfrl 

^ / H(t)dt = apKM / ^-^dr. (81) 

^0 Jo ifBj 

Solutions 

The procedure for solving this case is similar to the above case. The results are 
and 

T(r) = ' ( ] (m ^ = (83) 

where in terms of the parameters of the model y = [5n'^ / {{9r]'T')^^^][\^ / {aPK^^{Nr + 
2KMN)N'^g^)Y^^{LpBi/M)^ and in terms of the expansion e-folds from Eq. (^) 

From the specified conditions of temperature expressed through a, /3, and km two parameters 
in the model are determined 



^7r6 / f3N^{Nr + 2KMN) 



and 
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M ^'^ ^ \ A2 ) 3(r7T')V2 aV2(i - /3)Ari/2^ 

(86) 

Based on these solutions, some useful expressions are as follows. The number of mass sites 
that are thermally excited at a given instance are nt,e. = 2km/(|/t + l)^''^- The number of 
mass sites that ipo crosses during the warm inflation period is 

^ 3(Vr')V2 ai/2(i_;3)iVi/2 ' ^'"> 



with 

27r2 (3'>/^Ki^Ne{Nr + 2KMNY'^ 



™" 3(r/T')i/2 aV2(i _ /5)iVi/2 
The mass of the particles is 

2 Xvl 9iT'a{l-(3)KMNg^ 



(88) 



{yr + lf'". (89) 



For examining the adiabatic conditions 

d^o 3a(l-/?)«M V^o(r) vr^ a^/^(3"''4^{Nr + 2KMNfl^ 



47r2 AT^ 
The ratio of vacuum to radiation energy density is 

(r) 4 (3N, 



dr ANe iyr+l) 2{r]T'y/^ m/^g{yT + 1)^^ 

and 

rx(y) _ 

r^(M) " (|/r+ 1)2/5- 



T(r) (91) 



(92) 



Consistency Conditions 

The parametric constraints from the consistency conditions are the y^o adiabatic condition 
Eq. ([18[) with sufficiency condition 



Ne>-, (93) 
the thermal-adiabatic condition Eq. (|TUp 



TT^ / a/3(iV, + 2/s:MAr) \ ^ ^ 



2(r/T')i/2 V «:mA^ ; 
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(94) 



the force condition Eq. 



and the infrared condition Eq. 



AT' aNeig^ + p/sy 



> 1. (96) 



C. Discussion 

In this subsection, some general comments are given about the four cases examined in 
the previous two subsections. The first noteworthy observation is that for arbitrary e-folds, 
A^e; the consistency conditions impose very mild restrictions on the parameter space in all 
four cases. Although, the precise restrictions vary amongst the four general set of 

restrictions that is valid for all four cases can be given. First recall that by construction of 
the solution, the thermalization condition Eq ([T7|) always requires simply that a < 1. The 
most stringent restrictions arise from the force Eq. ( |2l]) and thermal-adiabatic Eq. (|19|) 
conditions, which together require 

< a ~ — . (97) 



(3kIi 100{Nr + 2KMN)' 

The other two consistency conditions impose very mild restrictions. The y^o-adiabatic con- 
dition Eq. ( p!8D is always accommodated provided Nf, > I and the infra-red condition Eq. 
(1^) is accommodated provided l3KMNg^/{aNe{g^ + P/S)"^ > 4 x 10"^. With all the con- 
ditions combined, they are fairly unrestrictive to the parameter space. As such, it leaves 
considerable freedom for treating density perturbations. 

Another interesting point is to compare the solution in Subsect. |1V B 1| for Tei ~ Tbi 



with those in Subsect. [IV A| for Tei ^ Tbi and see how well they match at some intermediate 



Tei < Tbi- There clearly should be some overlapping region since the model is exactly the 
same and only the treatment of temperature dependence is different. By comparing the 
expressions for A and ipBi/M, the functional dependence on the parameters is seen to be 
exactly the same except for with respect to (3. In regards to /3, both A and (pBi/M for the 
cases m Subsects. IIVAII , A^^ = 0, and |IVA2| , A^'^ > 2T/M, equate to the expressions m 



IV B 1| a.t (3 = 0.42 and 0.53 respectively. This is close to the approximate cut-off criteria we 
gave in 3B of /3 ~ 0.5 Eq. (0). 

The final point is that cross comparison amongst the four cases indicates several similari- 
ties amongst the solutions. In the remainder of this subsection, the origin of these similarities 
are examined. Alongside this, a qualitative understanding of the solutions are developed. 

The two basic equations of warm inflation, the (yjQ-equation of motion Eq. (^ and the 
energy conservation equation Eq. (0) have five properties, which are listed below. From 
the properties, all the similarities amongst the solutions then are explained. 
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property 1: The first time derivative of y^o, dipo{T)/dT, is related to some product of 
V9o(T)'^T(r)^ in both basic equations of warm inflation, Eqs. (P) and (0). (In the v^o-equation 
of motion, Eq. (j^), recall that we approximate the dissipative coefficient as I]*''^'(v?o ~ 
Mj)^?7(T) ^ T^ri{T)/{g'^M), where for the T^(t) term, in Subsect. A the time dependence 
is treated and in Subsect. B it is treated as a constant T^{t) = T^(0). In particular the 
dissipative coefficient in both subsections depends on some power of the temperature T(r).) 
property 2: In the y^o-equation of motion, the direction of dipQ{T)/dT is always opposite to 
the sign of ^Po{t). 

property 3: Since the yjo-equation of motion is first order, the solution requires one initial 
condition which then sets the overall scale for both ^poi'^) and T(r). In our approach, this 
scale is set by the condition T(0) = Tbi = hmM. As such ipsi is then a derived quantity, 
property 4: The force term in the (fo equation of motion in our approximation is always 

property 5: The number of e-folds is linearly related to the dimensionless time parameter 
at the end of warm inflation tei in all four cases oc ar^/. 

Properties one and two imply that the solutions have the general behavior 

cpoir) = ifBiDir) (98) 

where 

D{r) = (^:^^( or exp(-yr)) (99) 

and 

T(r) = f{\,g,f,N,^M,(y,P)y^rir)M'-^^ (100) 

with ■jfjjy'jT > and f{X,g,f,N,K,a,P) a function of the the model and thermodynamic 
parameters. Combining this deduction with property 3, we also can conclude that in general[| 

T{t) = f^MD-'^r) (101) 

The noteworthy point for the present discussion is the solutions for (Po{t) and T(r) are always 
a product of a mass-dimension one function which depends on the model and thermodynamic 
parameters and a time dependent decay function. The latter we represent through D{t) 
taken to some positive power with -D(O) = 1 and D{t > 0) < 1. For this discussion it is 
useful to think about the solutions this way, since the general features are contained in the 
mass- dimension one function. As such, the discussion to follow is not detailed about the 
decay function. 

The general behavior of the temperature permits two deductions. First from Eq.( |101D 
and the definition of (3 it follows that in general 



3t 



Equating Eqs. ( |10[)| ) and ( 101 ) lead to one of the two parametric constraints that in the previous 



two subsections determined A and lpbi/M. 
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y - (102) 

where h{P) depends on the specific time dependence of T(r) 0. Second, our approximation 
for the dissipative coefficient in all the cases studied has the general form 

E(^o - M,r^ ^ ^ ^ ^^D{rr-, (103) 

TV I TV MI TT 

where 7,, depends on the specific treatment of the temperature's time dependence. In Sub- 
sect. A 7.^ = 27t and in Subsect. B 7.^ = — 7t- 

By this point it should begin to appear evident that the differences in the various cases 
emerge primarily in the exponent of the decay function. This becomes fully clear once the 
basic equations are examined below. In fact it can be recognized that the y^o-eQuation of 
motion Eq.(P) actually is a equation of motion for D{t). Using Eq. (|103|) it has the general 
form 



On the other hand, from the general form of the ^Po{t) solution Eq. (|98|) , it also follows that 

oc (105) 

where = 1 if D{t) ~ exp{—yT) and zero otherwise. Based on this relation, we find 
that dp^^ipo) / dr = X(pQ{T)dipo{T) / dr /6 oc yp^{ipo)/{yT + 1)-'^^^'= which applied to the energy 
conservation equation Eq. (p!3|) implies 



r 



= ?MW ^ H(r) yr+l ^NM 

P,.(t)- 3t,^[N,. + 2kmND-'^(t)\T*(t) T^(M) y h(l}y"^' 

where jp/ and ■jp" are more exponents that depend on the specific solution regime. In 
Subsects. A and B 7p/ = 7^ and respectively. 7^" represents the time dependence of H. 
Thus 7p/' = in all the cases except the last, IIIB2, where 7p// = 2. 

Now we can deduce all the general features of the solutions in the previous two subsec- 
tions. In addition to examining the solutions in terms of the convenient parameters used in 
our analysis, it is interesting to see how the solutions depend on the parameters with direct 
physical interpretation, the Hubble parameter H and the slope parameter for (Po{t), y in 
exchange of a and N^. For p^{t)/ pr{T), the general expressions already have been given in 
Eq. ( 11061) . From Eqs. ([TO^ and (|10|) we find 



2T2(r) ^ ^ ^ ^ (yr + 1)1-^^ N, {yr+iy-^^' ^ ^ 



^Equating this expression for y with the one obtained with respect to the parameters of the model 
yields the second of two conditions that determine A and (pBi/M in the previous two subsections. 
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From Eqs. (|0|) and we find 

{H/r){Nr + 2kmN) ^ ' pNl{Nr + 2kmN) 



A ~ hirr , „ .Jf'' ■„ = Wrr t (108) 



and 

^ 1 ( {H/T)kI,{N, + 2kmN) \ ^ 1 /3V^4fiV,(iV, + 2/tMiV)V^ 

(109) 

Regarding the consistency conditions, the general forms are as follows: y^o-adiabatic condi- 
tion Eq. 

JL ^'^^M ^ ^M/?) D-"{r) ^ . 
the thermal-adiabatic condition Eq. (|TUp 

I, r]T'KljN ) {yr + iy-^^ [ {r]T')KMN J {yr + ly-^'^ ^ 

(111) 

and the force condition Eq. ( pT) ) 

(^r)V^ |^ (^/n(iV^+2/.MiV) y^^^3-3..(^) ^ j^ «/^^M(iV^+2/.MiV) y/^^3-3..(^) 

>0(1). (112) 

These expressions contain interesting insight into the solutions. The expression for 
Pv/Pr oc V{ipo)/T^, Eq. (|106|) , indicates that the relative content of radiation to vacuum 
energy is directly proportional to the Hubble expansion rate Hbi and inversely proportional 
to the slope parameter y. These general trends can be explained. The former is expected 
since for a slower the Hubble expansion rate, the red-shifting of radiation is slower. The lat- 
ter follows since as the slope parameter decreases, the decay of vacuum energy to radiation 
also becomes slower. 

The expression for m|/T^, Eq. ( |107| ), also has a simple explanation. It is independent of 
the Hubble parameter, which is indicative that this expression is entirely an outcome of the 
ipo equation of motion. Since the potential is a monomial in (po and the (po equation of motion 
is first order and overdamped, Eq. ( p.05| ) follows. Furthermore, in our case the ipo equation of 
motion always has two powers of the temperature in the form dipo/dr oc {dV {ipo) / dipo) /T"^ . 
Comparing this with dipo/dT ~ ypo from Eq. ( p.05|) and noting that our potential is quartic, 
V{p)q) ~ Ay^g, the expression for \p)'^/T'^ in Eq. ( p,07| ) follows. 

Based on these two expressions, Eqs. ( [L06|) and ( |107|) , the remaining expressions follow. 
For example since V{(Pq)/T'^ oc Hbi/v and {dV {p)q) / dp)Q) / {p)qT'^) oc y, it must follow that 
V9q/T^ oc HBi/y"^- It then implies that for a quartic potential A oc /Hbi- For a hypo- 
thetical monomial potential ^^(v^o) ~ XM'^~^(j)"', it still follows that p>q/T'^ oc Hsi/y^ but 
now X{p>o/My~'^ oc y^/HBi- The thermal adiabatic condition, oc {dp^/ dr) /T , also can be 
understood since {dip^/ dr) /T oc yp:>o/T oc H^bi ■ The behavior of the other expressions as 
well as other parameter dependencies can be deduced by similar reasoning. 



28 



V. ESTIMATES FOR AMPLITUDE OF DENSITY PERTURBATIONS 5p/p 



In this section estimates are made of scalar metric perturbations for flat spatial geom- 
etry, = 1, for the four cases in the previous section. For the calculations to follow, the 
basic formulas can be deduced either from first principles arguments or from a perturbative 
quantum field theory calculation, and we will take the former route. This approach has lim- 
itations, within which it is indisputable. For our purposes these limitations are unimportant 
in obtaining the basic formulas and in making the estimates for density perturbations in 
the quantum field theory model studied in this paper. The limitation of this approach is in 
sacrificing some of the details that could be obtained by the perturbative calculation. These 
details, which will be described at an appropriate place below, are not useful for our present 
purposes. Nevertheless, at some stage, it will be important to understand these details and 
to do the perturbative calculation. For this purpose, the present approach also is useful as 
a cross-check for the perturbative calculations and in providing a necessary outline of the 
formal problems that must be addressed for that derivation. 

To estimate density perturbations, the classical background field now is treated with 
fluctuations, (^(x, t) = ^poit) + (5v?(x, t), where ^poit) is the zero mode and 5v9(x, t) are small 
fluctuations about the zero mode. The calculations to follow are in momentum space. The 
Fourier transform of the fluctuations is defined as 

(5^(k,t)= / c/^x(5¥?(x, t) exp(-ik ■ x) (113) 
Jv 

where (k, x) can be either comoving coordinates or physical coordinates at a particular 
time. Hereafter we denote comoving coordinates as (k^, x^) and physical coordinates at 
time T as (kp(r), Xp(r)), where the r- dependence sometimes is not shown explicitly. The 
comoving coordinates can be regarded as the intrinsic labels for the modes of the scalar field 
fluctuations. Thus the modes will be denoted as (5(^kc ('?")• For definiteness, the comoving 
and physical coordinates will coincide at the end of warm inflation tei so that for any 
component direction ^^(t) = kc/R{T — Tei), where by our convention -R(O) = 1. For the 
four cases treated in the last section R{t) ^ exp[if(r)r]. Very often it will be necessary to 
consider modes in terms of their physical wavenumber at a given time r. Thus, the modes 
also may be denoted as 6ip{kp,T), where the (redundant) kc subscript has been dropped. 

Let us obtain the equation of motion for the fluctuations in flat nonexpanding spacetime. 
We argue below that this equation also suffices for our purposes in an expanding background. 
For the nonexpanding background case, no distinction is necessary between comoving and 
physical wavenumbers; the modes are denoted simply as 6ip(k,T). 

As stated earlier, the equations of motion for the fluctuations can be obtained by two 
approaches, either by general first principles arguments or a direct perturbative calculation, 
and here the former approach is taken. The basic fact that permits the former approach 
is that the calculation of the zero mode and fluctuation equations of motion, which follow 
the methods of p^-p^,p!^, are done in a near-thermal-equilibrium approximation, thus re- 



specting the fluctuation-dissipation theorem. This fact immediately determines the form of 
the fluctuation equation of motion, given that for the zero mode. In particular, assuming 
the fluctuations are small, their equation of motion is obtained through a linearization of 
the zero mode equation of motion, and with the inclusion of the gradient term, V^^v?, and 
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a noise term that represents the short distance dynamics of the heat bath. The fluctuation 
dissipation theorem uniquely determines the correlation statistics of the noise. The missing 
detail from this first principles procedure, which were eluded to earlier, is the relation of the 
noise function to the basic field variables of the Lagrangian. For our present purposes, this 
relation is not needed. Examples of such relations obtained from perturbative calculations 



are given in [|35| , p7 | 



Turning to the equation of motion for the fluctuations, as stated above, it is obtained 
through the linearized deviation of the zero mode equation of motion, with account for 
the gradient term and noise. Although our interest in the zero mode dynamics is in the 
overdamped regime, where its equation of motion becomes first order in time, for obtaining 
the fluctuation equation of motion we must start with the initial second order zero mode 
equation. Compared to the zero mode equation of motion, the only additional term in the 
fluctuation equation of motion that could make it second order is the gradient term. This 
will happen if the gradient term creates a sufficiently large curvature (= — + m'^(j)) in 
the potential to overcome the large damping force term. However, this will not happen 
since our interest is in fluctuations that relatively are not very large. In particular, we will 
see below from explicit calculations that for the fluctuation wavenumbers of interest, the 



curvature term (= y + m^) generally is smaller than the dissipative coefficient (recall 
this is the term multiplying dipo/dt in Eq. (|^)), and generically for a damped harmonic 
oscillator equation of motion this implies the overdamped (first order) regime. Thus, the 
equation of motion for the fluctuation 6ip{\<., r) is obtained from the equation of motion for 
(Po{t) that is given in Subsect. [Ill A| Eq. (^ and derived in [p!6| . 



T(y,o, ^)^^^^ = -(k' + mli^, T))S^iK r) + ^(k, r) (114) 



where 

t.e. 



T(^o,T) = ^(^0 - M,fr^,{T)V^{M) (115) 

i 

and m^{{pQ,T) is given in Eq. (^. This equation is obtained from Eq. (^ by substituting 
<y9(x, t) and retaining terms linear in Sip{'K, t) and Fourier transforming to k-space. In addition 
a noise function ^(kp,t) is added which respects the fluctuation-dissipation theorem 



<e(k,t)>5=0 (116) 



< e(k, r)e(-k', r') >/=°° 2T(^o(r), T(r))r(r)(27r)3<5(3)(k - k)<5(r - r'). (117) 

The equation of motion ( |1 14| ) represents the standard near-thermal-equilibrium dynamics in 
which ^(k, r) drives the correlations of 5v9(k, r) to thermal equilibrium with the relaxation 
rate of the initial conditions (k^ + rn?^{(.pQ,T)) /T {(.p ^T) . 

The relevance of Eq. ( [L14| ) to expanding background is that it approximates the equa- 
tion of motion for a mode 5v^kc('7") during a Hubble time interval say Xj ~ r ~ Tj+i 
with Tj+i — Ti r^{M)/H. Furthermore, the momentum vector in the equation of mo- 
tion Eq. ( |1 14| ) is identifled with the physical momentum of k^ which to a good approx- 
imation is flxed at one intermediate time during the respective time interval, such as r,, 
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k — s> kp = kcexp[— i7(r)(rj — tei\- This approximation is valid for large kp, when the evo- 
lution of if(r),kp(r) and ^{j) is adiabatic relative to the evolution of 5ip\f^X^) during the 
respective time interval. Within the regime where this approximation is valid, the evolution 
of 5ip\f^X^) can be computed through piecewise construction of solutions for a sequence of 
Hubble time intervals, similar to the demonstration in |]1T|. We will see below that the 
regime of kp where the above approximation holds also is the appropriate regime for our 
purposes of estimating density perturbations. 

Consider what happens to a mode i^v^kcl'^) that is immersed in a heat bath and is in an 
expanding background. The larger k^ is, the faster is the relaxation rate. If k^ is sufficiently 
large for the mode to relax within a Hubble time, then that mode thermalizes. Thus at any 
instant during expansion, one can expect modes with physical momenta bigger than some 
lower bound kp to thermalize within a Hubble time interval. For these modes, within a 
single Hubble time interval, the flat-space equation of motion for the fluctuations Eq. (|114|) 
is approximately valid. 

As soon as the physical wavenumber of a v?(x, r) field mode becomes less than kp, it 
essentially feels no effect of the thermal noise .^(kp, r) during a Hubble time. Thus for mode 



S^kci^), it essentially does not change once 



|kp| — 



|k,|exp[-(if(r)/r^(M))(r-r^,)] < 



kp, and at |kp| = kp the mode assumes its thermalized distribution. If kp > H{t), it 
implies the Sip^^ (r) correlations that must be computed at time of Hubble radius crossing, 
|kp(T)| = H{t), are the thermalized correlations that were fixed at |kp(r)| = kp. This 
effect was clarified for warm infiation by Yokoyama and Linde |^5| and they referred to it 
as "freeze-out". 

In order to determine kp, consider the solution of Eq. ( |114| ) for 5(y9k^(r) within the 
Hubble time interval tq < t < tq + T/H{tq). We will ignore the time variation of kp, 
T{ipQ,T), and m^{ipQ,T) during this time interval. Their values at tq will be used kp(ro), 
T(ro) = T{(po{TQ,T{ro)), and ^^(ro) = m^{Lpo{TQ),T{To)). The approximate solution then 
is 

S'^k^ir) ^ — — exp[ ^ t (r - tq)] / exp[ 



T(ro) 



k^ 



T(ro) 



TO 



T(ro) 



T(ro) 



(r-ro)]. 



romk„r')dT' 

(118) 



and for the corresponding correlation function 



1.0 — exp 



T(ro) 



{r - ro) 



+ ('5</'kc(^o)'5</'k^(^o))?exp 



2(kg + m|(ro)) 
T(ro) 



{r - To) 



;ii9) 



When the exponentially decaying terms are negligible, the above correlation is equivalent to 
the high-temperature correlation function {5ip{kp,TQ)5ip{k'p,TQ))T- 

In this solution Eq. ( 118 ), on the right hand side, the first term is the noise contribution 



that is driving the mode to thermal equilibrium and the second term contains the memory 
of the initial conditions at r = tq, which are exponentially damping. By definition of freeze- 



out, for |kp| ~ kp the second term damps away within a Hubble time and for |kp| ~ kp it 
does not. To quantify the criteria, the freeze-out momentum kp is defined by the condition 
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k| + m^(y.(r),r(r))r,(M) ^ 

T(^o(r),T(r)) i7(r) " ^ ^ 



This relation allows us to follow-up our discussion above Eq. (|114|) about the fluctuation 



equation of motion being first order in time. From the above equation we see that the 



curvature term in the potential + is less than the dissipative coefficient (in this 
notation comparing to Eq. (^ the dissipative coefficient is T/T^) since generically H < 
T/r^. As such, the overdamped equation of motion for the fluctuations, Eq. ( |114| ), is 
justified for the wavenumbers of interest to us. 

Now we can write down the basic equations for calculating density perturbations during 
warm inflation. Properly this should be fully derived from the linearized General Relativity 
equations for perturbations [^8| -|50[]. Nevertheless, an examination of the adiabatic density 
perturbations derivation for supercooled inflation [^,0 with the modification of a subdom- 



inant radiation component for the warm inflation case leaves unaltered the basic expression 



for 6p/p. This conclusion was arrived at in our earlier papers |[T^JTO[| and independently by 
Yokoyama and Linde P3[. Therefore from this we find 



S_P(. ^ V'{^o)AvHir)[kF{r-f{T))] ^ 6H{r)AvHir)[kF{r-f{r))] 

^ {T^{M)d^o/dTy + {4/3)pr{T)^ 5T^{M){d^o/dT) ' ^ ^ 



The middle expression is the one used by Berera and Fang |jT2|| and the latter is the 



Guth and Pi expression |^D[. During warm inflation, since from Eq. (|T3]) pr ~ 
V'{ipQ)V^{M){dipo/dT)/{AH{T)) and pr > {V^{M)dipQ/ drf , the middle and left expres- 
sions above are equivalent up to an 0(1) constant. It would be interesting to in- 
vestigate which of the two expressions is fundamentally more proper. In Eq. (|121D , 
Ai{)H(T)[kF{.T — ^(t)] is the amplitude of the scalar field fluctuations and it is composed 
of all the modes whose physical wavelengths cross the Hubble radius within one Hub- 
ble time interval about time r. Recall that the comoving mode of physical wavenumber 
|kp(r)| ~ H{t) had its amplitude frozen at an earlier time r' when its physical wavenumber 
was |kp(r')| ~ A;^(r'). Since |kp(r')|/|kp(r)| = exp[(/7(r)r - /7(r')r')/r^(M)], this im- 
plies kpir') = H{T)exp[{H{T)T - hIt')t')/T^{M)]. This defines f(r) = t-t'. For the 
cases in the previous sections, kpir) is slowly varying, so we will use the approximation 
kp{T — fij)) ~ kpi^r). The expression for the scalar field amplitude is defined as the natural 
finite-temperature extension of the T = expression of supercooled scalar field inflation and 
with account for kp- We use the definition 

A^llkp] ^ —fj ciS(5^(x,r)5^(0,r))^exp(-^k, ■ x,) (122) 

where the kp - shell is defined as the spherical shell which is bounded between kpe~^f'^ < 
|kp| < kpe^^'^ (we approximate the shell thickness simply to be kp). The expression on the 
far right in the above equation is valid when kp < T. The definition Eq. ( |122| ) is equivalent 



to the one given by Linde ||5T]] in which one retains from < (x, r) >t the contribution 



from wavenumbers within the kp-sheW. From Eq.(7.3.2) and Eq. (3.1.7) of Linde's book ||51 
this gives 
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(123) 



When kp > T, Eq. ( |122| ) implies Aipjj{kF)kF>T = kp/i.'^T^ )■ For this region, it is a 
poor approximation to use the zero-mode dissipative coefficient in the k-mode equation of 
motion Eq. ( |1 141 ) . For this regime, a proper calculation of the k-mode dissipative function is 
important. Our calculations in the next two subsections consider only the high temperature 
expression on the right hand side of Eq. ( 122 ). 

Substituting the expression for /S.ip'\j from Eq. (|122|) into Eq. ( |121| ), the final expression 
for 5p/ p for the mode that crosses the Hubble radius at time r during warm inflation ig0 



5p u,<T m{r)/T^{M)]^kF{T)T{T) 
— (r) ^ -= 

P h\/2'KdipQ{T)/dT 



(124) 



IS 



Recall, the comoving mode that crosses the Hubble radius at time r 
H{T)eM{H{r)/T^{M)){T - TEi)]. 

From the above considerations, the final prescription for computing 5p/ p is simple. First 
determine the freeze-out wavenumber from Eq. ( [1201) , and then substitute this in Eq. (|124|) 
along with expressions for T(r), H{t) and dipo{T)/dT from the previous section. In the next 
two subsections this is done for the cases in the coinciding subsections of the last section. 
Finally Subsect. |V C] follows with a discussion of the results. 



A. Tbi > Tei 



This subsection considers the cases in Subsect. A of the last section. 



1. Nr = 



From Eq. (|114|) and taking T in Eq. ( p.l5| ) the same as in Eq. 
motion for the fluctuation is 



d^V9(kp, r) 
dr 



-[^ + mJ(v.,T)]-— 



TT 



-5v5(kp,r) = -2ly 



\ Vbi ) 



2/7 



the equation of 



k2 

^ 5v9(kp,r). 



(125) 



^In our estimates of density perturbations, the equation of motion for 5ip\f_^{T), Eq. ( [L14D uses the 



k = dissipative function Eq. (|lO|). The derivation in |37,1^ suggests that the dissipative function 
?7i(T) decreases as k increases, since it requires "off-diagonal' Green's functions. This would imply 
kp decreases since the relaxation rate is faster, thus any mode (^(/'^^(t) would thermalize faster 
at every physical wavenumber. Ultimately this means for a given kc mode that is crossing the 
Hubble radius, 5p/ p decreases relative to our estimates. In the future, a proper derivation of the 
dissipative function at non-zero k-vector would be useful. 
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The expression on the extreme right was obtained by substituting the expression for y below 
Eq. ( p8l) and considering the regime where ^ rn?^{ipQ,T) ^ \lPq{t)I2. The freeze-out 
wavenumber from Eq. ( |120| ) is 

Using this, from Eq.( 124 ) we find 

3/4 



P 



a^"(3^'^K]!,'N'l'g^l\yT + 1)^/^ (127) 



57r7/2 

The spectrum is very flat. Since the relation between time r and the comoving wavenum- 
ber kc crossing the Hubble radius at time r is r = — (r^(Af)/i7(r)) ln(|kc|/-ff£;/) + tei ~ 
-{T^{M)/Hbi) H\K\/HBi)+rEi, the variation oi5p/p cx [y{-{T^{M) / Hbi) \n{\K\/HBi) + 
Tei) + 1]^''^ is only logarithmic. For example, for arbitrary e-folds iVg, 5p/ p for the last scale 
that crosses the Hubble radius (smallest scale) is a factor bigger than the first scale 

that crosses the Hubble radius (largest scale). Observe that this deviation from exact scale 
invariance is an outcome of the nontrivial thermodynamics and is not, in particular, initially 
inputed by hand by choosing a nonanalytic potential. 

2. Nr > 2NT/Tei 



From Eq. (|114|) and taking T in Eq. ( p.l5| ) the same as in Eq. (24), the equation of 
motion for the fluctuation is 

^ _m + ,„J(,,r)l-,^,^^.,(k„.) ^ -6, (128) 

The expression on the extreme right was obtained by substituting the expression for y below 
Eq. ( ^51) and considering the regime where k^ ^ m'^{ipQ,T) ^ Xipl{T)/2. The freeze-out 



wavenumber from Eq. ( |120| ) is 



.^,^^^_{Hbi/T^{M))\vI{t) 



kUr) = " (129) 



Using this and from Eg. (|124|) we find 



p 7r7/2 0^ - 7rV2 ^1/2 • ^^'^"^ 

For this case, the spectrum is exactly flat, i.e. 5p/ p is the same for all e-folds. 

B. Tbi ~ Tei 

This subsection considers the cases in Subsect. B of the last section. 
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1. Vb > 

From Eq. ( |114| ) and taking T in Eq. ( |115D , the same as in Eq. (|5^, the equation of 
motion for the fluctuation is 

(131) 

The expression on the extreme right was obtained by substituting the expression for y below 
Eq. ( |65D and considering the regime where ^ m'^{ipQ,T) ^ \ipl{T)/2. The freeze-out 
wavenumber from Eq. (|120|) is 



4(r) = ^Wr)/rAM))X^Ur) ^ ^^^^^ 



Using this and from Eq. ( |124| ) we find 

Sp. . Q^{r]'T'f'^{HBilT^{M)fl^N^l^g^l'' 



57r7/2 (AT,, + 2KAfA^)V2 



{yr + 1)^/^ (133) 



The spectrum is nearly flat with deviations that are logarithmic. For arbitrary e- folds Ne, 
5p/ p at the last e-fold (smallest scale) is only a factor /^"^/^ bigger than at the first e-fold 
(largest scale). Recall that the regime for this approximation requires 1 > /5 ~ 0.5. 



2. Vb = 

From Eq. (|114|) and taking T in Eq. ( p.l5| ) the same as in Eq. (^), the equation of 
motion for the fluctuation is 

d5^{K^r) vr^T(r) ^ 9 ( Mr)^'" ^> , . 

dr r]T'Ng^KljM'^ 5 \ ^bi J ^"PoiV 

(134) 

The expression on the extreme right was obtained by substituting the expression for y below 
Eq. ( p3D and considering the regime where k^ ^ rn?^{ipQ,T) ^ X({)1{t)/2. The freeze-out 
wavenumber from Eq. (|120|) is 

miiEim^i^r^ ,135) 

In this case the Hubble parameter from Eq. (^) is r-dependent. Using this expression and 
the above condition for kp-, from Eq. (p.24|) we find 
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5n ^ QV2{iT'f" {Hsi/T^{M)f/'N^I'g^/^ 

~ {N, + 2kmNY/^ + • ^^^^^ 

The spectrum is nearly flat with deviations that are logarithmic. For arbitrary e- folds N^, 
5p/ p at the last e-fold (smallest scale) is only a factor bigger than at the first e-fold 

(largest scale). Again recall that the regime for this approximation requires 1 > /? ~ 0.5. 

C. Discussion 

Some aspects of the results for density perturbations are discussed here. Recall from 
observation that 5p/p ~ 10"^ From our results, for all four cases we find 5p/ p ^ 



10^'^a^^'^K%'^N^^'^g^^^ . Thus the parameters a, km, N and g must decrease 6p/p by at least 
three orders of magnitude. For observational consistency it is permissible if Sp/p < 10~^ 
from inflation, since post-inflationary mechanisms such as cosmic strings also can produce 
density perturbations. 



Overall consistency with respect to quantum field theory from Subsect. [1IIB| and obser- 
vation, Sp/p ~ 10~^ and Ne ~ 60, can be achieved in a variety of ways. From the discussion 
Subsect. |1V C| , Eq. (^) is a general consistency regime for quantum field theory for arbi- 
trary A^^e. Combining this with the general form for Sp/p stated above, overall consistency 
can be achieved for all four solution regimes in general for 

Alternatively, for arbitrary a, (3, and km, Sp/p can be made arbitrarily small by requiring 
g ^ 0, Ng'^ = const, with a < 6. Thus the model provides sufficient freedom to achieve 
consistency independently with respect to quantum field theory and observation. Further- 
more, both alternatives for decreasing Sp/p, a — or (7 — improve the validity of the 
underlying approximations. In the former case the adiabatic regime is deepened and in the 
latter case perturbation theory is further justified. 

It is important to note the magnitude of the Hubble parameter within the consistency 
regime. From Eq.(||) recall that Hbi = a/3KMT^{M) ^ [a/?KM/(1927r)](c/^ + /VS), with 
thermalization condition Eq. (|1^) requiring a < 1. None of the consistency conditions 
so far have imposed any conditions on the self-coupling parameter / except perturbation 
theory which requires / ~ 1. Thus irrespective of g, the bound can be given Hbi < 
(2 X 10-^)a/3KMM. 

There is some freedom to this bound. Observe that the entire calculation is independent 
of the number of bosonic decay channels for the respective Xifc-field. This refiects itself in the 
coupling constant dependent factor C in r^(T). Note that in the basic equations the product 
?7i(T)r^(M) always arises together, and because of this the coupling constant dependent 
factor, in our model C ~ {g^ + /^/8), cancel. This is not a coincidence, since fundamentally 



rji (T) ~ l/F as evident from the formalism of [p7|,p^ and from heuristic arguments in 
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For the model in this paper, each Xik field has only two decay channels, one to itself and 
the other to the field. The Xik field also may interact with other fields. For example 
consider the additional interaction ((7^/2) X^j^i X?fcC"|; where aj are bosonic fields that only 
interact with Xik- In our notation, these cxj fields are nondissipative heat bath fields^. Such 
interactions modify the Xik decay width as 

-p ,rr\ add a interactions T ^ f ,^ 4s /1 ooN 

Tx.^u) ^ 192^^^ +^ + Nog„). (138) 

If all the couplings (7^ > 0, then one must restrict the number of a fields, N^, since their 
interaction with Xik also adds a thermal mass contribution g^N^T'^ /12. In order to keep 
m^^^. < T, it requires g^N^ ~ 12. In the strong coupling limit g^ ~ 1, this modification 
with its stated limits increases r^(M) by a factor 10^. In this case the bound on the Hubble 
parameter is Hbi < (2 x 10~'^)al3KMM. 

To go one step further, the thermodynamics of warm infiation in the previous sections 
would be unmodified if in total as many as N^r ~ i3kmN a-fields coupled to all the Xik fields 
that are thermally excited at a given instant. In this case Hbi < (2 x 10~^)a/5^fi;|^A^M. For 
this case, precaution is necessary to keep m^-^, < T. One way to achieve this is if the Xi ~ 
couplings are not sign definite. Since the one-loop thermal mass correction is sensitive to 
the sign oc i^f^ whereas oc g^ is not, in principle F^ can be made arbitrarily large while 
the growth of the thermal mass correction is controlled. These modifications to F^ may be 
useful for phenomenological applications. 

The final point addressed in this discussion is the general behavior in all four cases of 
the density perturbation formulas. The equation of motion for mode 6ipk^(kp,T) can be 
obtained from a linearized approximation to Eq. ( |104[ ) and by replacing A(/?o which 
gives 



ddifCkp, t) k. 



dr r]T'Ng^Kl 



-D^^{T)6^{kp,T) + noise. (139) 



M 



The decay function -D(r) = {lpo{'t)/vbi) and all exponents 7 are the same as defined in 
Subsect. |IV C| . Once again, our focus is not on the (slowly varying) time dependence. 



but we have included the correct r-dependence for completeness. The quantity multiplying 
6ip{kp, t) on the right-hand-side of Eq. ( |139| ) is the decay rate for the respective mode. Thus 
the freeze-out momentum for the modes is 



Hir) V" 



{7]'T'fl^N^'^gKMD^^-+^p"^'\T)M 



~ {iV'f'^a^/^P^'^4^N^/^gKMD^^-+^^"^/\T)M. (140) 

Using this expression and the general expressions for dipQ{T)/dT, Eq. (|105| ) and (psi/M, 
Eq.O, implies 



Such secondary interactions with (p also induce dissipative effects, though we will not treat that 
here 
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Some features are worth mentioning. In Eq. (|140|) , kj^ should be proportional to the 
Hubble expansion rate Hbi, since a slower expansion rate allows longer relaxation time, 
thus modes of lower |kp| can equilibrate. The dependence of Sp/p on Hbi in Eq- ( |141|) is 



less than the naive linear behavior given by the defining formula Eq. ( |124| ). This arises 
because of Hbi dependence induced by the dynamics on the other factors, oc H]^^^ and 

(T/cpBiy/^ oc H^y\ 



VI. EXAMPLES 

In this section the solutions from the last two sections are studied in a few examples. 
It is not the purpose of this paper to detail the phenomenological consequences of these 
solutions. However, here we would like to obtain some idea about the absolute scales for 
the various dimensional quantities and how they chance in various parametric regimes as 
well as in the limit of increasing adiabaticity. In much of this section, the numerical value of 
premultiplying constants are evaluated in various expressions and we set i]' = 48, T' = 1/192. 

Before turning to the examples, note that in our construction, since the Hubble pa- 
rameter, Eq. (|T^) , is proportional to r^{M), the scales of all warm inflation quantities are 
controlled by this decay rate. Recalling our comments from the discussion Subsect. [V Q 
about modifying C in r^(M), in "our model" C = + /^/8 and in the "extended model" 
considered in Subsect. |V C| where a-fields were introduced C = + p/8 + g^N^j. In the 
three subsections that follow, we quote estimates for both "our model" and the "extended 
model" . All the estimates to follow always are within the observationally consistent regime 
with respect to expansion e-folds N^. ~ 60 and density perturbation Sp/p ~ 10~^. 



In Subsect. |VI A| a general estimate of the scales M,Tbi, HBijTrifpgj, and iminM are 



given that encompasses the four cases studied in the last two sections. Subsects. [VI B| and 



VI C] focuses on the case studied in Subsects. [IV B 2| and [V B 2[ for Vq = and Tei ~ Tbi 



In Subsect. [VI B , the dependence of the scales and parameters are examined in the limit 



of arbitrary adiabaticity a — 0. Finally in Subsect |V1 (J[ the warm inflation solutions are 
examined as the infiaton self-coupling parameter A varies over a wide range including A ~ 1. 
Within the limits of the present analysis, we will find that A is restricted to be tiny. However, 
a possibility is examined that could increase A up to 10~^, within a regime that is consistent 
with observational requirements on e-folds and density perturbations. 



A. Estimate of Scale 

The absolute scale of all dimensional warm inflation quantities in our solution are deter- 
mined once the Planck mass is introduced Mp ^ lO^^GeV. Our solutions are constructed 
such that all the dimensional quantities have been expressed in terms of M. To determine M, 
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n ote that the H ubbl e parameter Hb i can be expressed in two way, by its definition Hbi = 
^87rp(0)/(3m2) ^ ^87rp„(0)/(3m2) and the expression Eq. Hbi = af3KMT^{M). For 
Pt,(0), from the expressions for the ratio Pr/pv Eqs. (]36D,(|5^), (|72[), and (pO|), note that 
p^(0) ^ A/'ePr(0) = NeTi'^iNr + 2kmN)T^j/16 = Ne7i'^{Nr + 2km N) / 16. Using this in 
the first expression for Hbi and Eq. (|I^) in the second and equating the two, we find 

M = , ^ ^Mp. 142 



An approximate upper bound can be obtained using Eqs. (^T]), (|57|), ([f7|), and (Q) which 
imply a ^ /3/[200(l + Nr /{2KMN))] and Eqs. (|2p, (0), and which imply 

Km ~ 10-^1 + Nr/ {2kmN) / (3. Substituting these bounds into the above equation Ol42|) 
gives 

(8 2X10-^),.V^C 
y]V3V(l + iV,,/(2K„JV))9/' ' ' ' 

From this expression and the bound on km it follows that 



(8.2 X 10-^)/?5/2(- 



< (6.8 X 10-^2)/39/2C^ 

"""^ ^ VAyV(l + iV./(2.MiV))-/4^- (145) 

< (4.4xlO-^)/3V2(l-/?)VV 

~ iV.(l + iV./(2.MiV))^ "^-^ 

and g^Ei or equivalently the scale of the mass levels ~ Vin^^ in the DM-model is 



;i.o X io-^)/?=^±^^"iv,i/2^ 

(1 - /30(i))Ari/2(l + Nr/{2KMN)fl^' 



To quote some numbers, consider a typical case with (3 = 0.5, A^'g = 60, Nr = 0, 
= 5 and "our model" with / ^ 1 so that C ~ 1/8 ("extended model" with g^r ~ 1, 
A",, ^ 12, so that C ^ 12). For this case we find M ~ 5.2 x 10®GeV (~ 5.0 x lO^^GeV), 
Hbi ~ 2.7 x lO^GeV (~ 2.5 x lO^GeV), Tbi ~ 1-0 x lO^^GeV (~ 1.0 x lO^^GeV), z^inM ~ 
1.2 X lO^^GeV (~ 1.0 X lO^^GeV), and Tei = Tbi/2. We also find A^^ ~ 2.6 x 10^ mass 
sites are crossed. If we set km at its lower bound, km = l/y/2aP, and require 5p/p < 10~^, 
then in all four cases within the high temperature regime, kp < T, it requires g < 0.2. From 
this it follows that A < 10~^^ and m^gj ~ (2.1 x 10~^)Tb/. The high temperature validity 
regime for the density perturbation results in Sect. [V| require kp < T. From the expressions 
for kp in the four cases, we find in general kp ~ y/W^m^. Thus the estimates given here 
involving density perturbations are valid for < 2 x 10^. Finally, for both our and the 
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extended models, the thermalization rate r^(T) is about 400 times faster than the Hubble 
expansion rate H{t) ~ Hbi, so that the thermalization approximation is well satisfied. 

Although A is tiny, the infiaton mass, is large relative to the Hubble parameter. In the 
above example is three orders of magnitude below the temperature scale but four orders 
of magnitude larger than the Hubble parameter. The smallness of A preempts questions 
about fine tuned potentials, similar to the situation in supercooled dynamics. This point 
briefiy is addressed in Subsect. [VI Q . However, it should be noted that for these tiny values 
of A, when the thermal damping is removed after the mass site M;^.^, the potential does 
not support infiation. Once the thermal damping is removed, the only damping term that 
remains is due to the coupling of the infiaton to the background cosmology. This yields a 
3i?0o term that is familiar from supercooled infiationary dynamics. The infiaton equation of 
motion then becomes SHdipo/dt = —Xip^/Q = — m|(y9o/3. Thus in a Hubble time At = 1/H, 
\A{Pq/{Pq\ ^ m'^/{9H^) ^ 1, so that (po rapidly falls down the potential. In words, the 
curvature of the potential is huge relative to the scale of the Hubble expansion rate. As 
such, to terminate warm infiation in this model and go into a radiation dominated regime, 
it suffices simply to stop coupling the infiaton to mass sites. 



B. The Limit of Arbitrary Adiabaticity 

The self consistency of the near-thermal-equilibrium quantum field theory formalism ap- 
plied in this paper requires satisfying the conditions in Subsect. [Ill Ij| . Based on the solutions 
in the previous two sections, we find that the most constraining consistency conditions are 
the thermalization and thermal-adiabatic conditions Eqs. ( jlT]) and Eq. (|19|) respectively. 
As it turns out, both these conditions are controlled by a single parameter in our solutions, 
a, with the validity for both conditions improving as a — 0. In this subsection the obser- 
vationally consistent regime with respect to and Sp/p is studied as a function of a, in 
particular, in the limit of arbitrary adiabaticity a — > 0. We focus on the case in Subsects. 
rVB2[ and [VBl for Vq = 0, Tei ~ Tbi- 



For the case of interest, from Subsect. [IVB 2[ the consistency conditions are 



a < (148) 

and 

where throughout this subsection we set A^^ = 0. The mass splitting scale parameter M is 
determined by the same procedure as in the previous subsection. We find 

^ = TTTJ^ ./o Mp. (150) 

The other dimensional quantities can be determined easily from this. 

In Fig. 1, the a dependence of all dimensional scales are shown for two cases, km is 
set to its lower bound, km = 

l/(2y2aV2^i/2)^ ^j^h always the restriction km > 1- The 
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limit of arbitrarily increasing adiabaticity is a — > (— log^gl'^) ~^ ^^e scales, 

M,Tbi, HsiyirainM, and m^gj, are in GeV with their log^^Q plotted. The sohd lines are for 
the case iV = 5, A^'e = 65, /3 = 0.5 and C = 1/8 ("our model"). The dashed lines are for the 
case A^^^ = 1/8, = 65, /? = 0.5 and C = 1/8 ("our model"). For the "extended model", 
C = 12, all scales in Fig.l are shifted up by a factor ~ 10^ except for Hbi which is shifted 
up by a factor ~ 10^. For the region to the left of the dotted vertical line at 0.6 {a > 0.25), 
Km = 1. 

M, Tbi, Hbi, and Vin^ are independent of 6p/ p, whereas oc g, so it depends on 6p/ p 
since from Eq. ( p.36| ) 

a ^ 30 45-Mp}!!1_ (151) 

For a —>■ with everything else fixed, g increases. Since g < 1 is required by perturbation 
theory, the model requires Sp/p as a — > 0. In Fig. 1, we set Sp/p = 10~^ down to 
the smallest a possible, which is given by the vertical solid and dashed lines for the two 
respective cases. To the right of these lines (smaller a), 6p/p is less than 10~^. For the 
other parameters in the model, as a ranges from to 1 (— log]^Q(a) ranges from oo to 0), for 
the solid case the ranges are g from 1 to 0.013, A from (oc a) to 4 x 10"^'', km from oo 
(oc 1/a^/^) to 1, and Nm from oo (oc 1/a) to 350; for the dashed case the ranges are g from 
0.59 to 1 X 10-^ N from 1 to 3 x 10^^ A from (oc a) to 4 x 10~^\ km from oo (oc 
to 1, and Nm from oo (oc 1/a) to 350. 

For a < which is in the region to the right of the dot-dashed vertical line at 

— logio(a) = 1-9, all adiabaticity conditions are valid. To the left of this vertical line, the 
thermal-adiabatic condition in its stringent form is not valid. However, as discussed in 
Subsect. |111 C|, the thermal-adiabatic condition may still hold in some part of this region. 



To determine the extent to which the thermal-adiabatic condition can be relaxed requires 
details about thermalization that go beyond the simple high-temperature approximations 
applied in this paper. 

The DM-model warm inflation calculation in EHl is similar to the case in Subsect. liV B 2 



for Vo = 0, Tei ~ Tbi- The difference is ignores the minor modifications that arise due 
to time dependence of the temperature, whereas this was treated in Subsect. |1V B 2|. The 



region studied in is for 0.5 < a < 1. This was chosen for its simplicity in illustrating 
the basic features of the results. In this region, all the consistency conditions are satisfied 
except the stringent form of the thermal-adiabatic condition. Based on earlier discussions 
in this paper, this region is still within the plausible validity region. We have verified that 



in the region of overlap the results in |26| agree with thoses in Subsect. [VB2 



C. A Dependence of the Solution 

In this subsection, we examine the A dependence of our solution for the case from Sub- 
sects. |1V B 2| and |V B 2| . For this, we treat A as an independent variable in exchange for km, 
which from Eq. (|85|) gives 

/53 t\t3\ 
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where a = (2/57r^)~^ has been set to its upper bound and we only consider the regime with 
Nj. = 0. Here and throughout this subsection, the numerical values of all constants are 
quoted and we have used T' = 1/192, r]' = 48. The parameters on the right-hand-side of 
Eq. ( |152| ) can be varied freely up to the mild constraints on km, Eq. (|149|) . Substituting 
Eq. (|152D into our solutions in Subsects. ( [IV B 2| ) and ( |V B 2| ) we find 



Tsj = (2.1 X lO-^ r M„ (154) 



Hsr = (1.8 X 10-^^)L_^|^M„ (155) 



m^,, = (3.9 X 10-^) ^ 1^^^ M„ (156) 



and Nm = (3 X 10^^)P^{1 - P^/^)XN^/[{1 - p^Ng^] so that 



^_3, (1 - Pf/Hl-P' /')g'C 



i^^M = (5.5 X 10-«)^ oMU^T ^ 12 (157) 



Also we express g in terms of A and Sp/p to obtain 

This can be substituted in the above expressions for the scales, but the perturbative restric- 
tion must be respected g < 1. In the observationally consistent regime, 5p/p ~ 10-^ this 
perturbative restriction is comfortably satisfied for a wide range of the parameters. 

It is interesting to examine the maximum size of A within the observationally consistent 
regime of warm inflation. For satisfying just the horizon/flatness problems, the A ~ 1 
regime has solutions if also (? ~ 1. However the primary restrictions arise from the density 
perturbation constraints. If we restrict to only the high temperature regime for the freeze- 
out momentum, the inequality kp < T must be imposed to Eq. (|135|) . Reordering this 



inequality to isolate A and using the solutions from Subsect. |IV B 2| , we find the constraint 

A<(6.6xl0-)ill^. (159) 

Setting A at its upper limit and substituting into Eq. (|158|) , for Nf. = 65, Sp/p = 10~^, 
p = 0.5, we find ^ ^ 7.2 x lO'^, which from Eq. ([T59| ) implies A ^ 1 x 10"^^ This 
value of A is two orders of magnitude larger than the limit in Subsect. [VI A| , because here 



Km ~ (1-4 X 10^) /N is not at its lower bound. With these values for A,^', iVg, /? , and leaving 
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unspecified, from Eqs. ( |153|) - ( |157| ) for "our model" with C = 1/8 ("extended model" 
with C = 10) we find the scales in GeV M ^ 0.7A^(6A^), T^/ ^ 1 x 10^ (8 x 10*^), Hbi ^ 11 
(7 X 10^), m^^, ^ 2 X 10^ (2 X 10^), and Nm = W^N so that VinM ^ 2 x 10" (5 x 10^=^). 

For kp > T, as stated earlier, the high-temperature calculations in the previous section 
are not valid. Suppose that thermal dissipation is inactive for wavenumbers larger than 
the temperature kp > T. In this case, if the freeze-out wavenumber from Eq. (|135|) is 
larger than T, the 93-amplitude should be set to its limiting value Ayj^^ k, T^/(27r^). In 
this case, the density perturbation expression in Subsect. [VB2| is replaced by 5p/p ~ 
3a^/^/?^/^(7(l/r + l)^/^°/(57r'^). From the previous paragraph, kp > T corresponds to the 
region of A in Eq. ( |159| ) with < replaced by >. 

For this hypothetical case, if we set A^^e = 65 and /3 = 0.5, we find ^ (1 x 10^^)A/A^, 
M ^ (8 X 10-^^)iVCMp/A3/2, Tbi ^ (9 x lQ-^^)CMp/ X^'^ , Hbi ~ (8 x lQ-^*)C'^Mp/\^/^, 
m^3, ^ (2 X lQ-^^)CMp, and Nm ^ (8 x W^)\/N so that Vi^M ^ (6 x lQ-^^)CMp/ X^'"^ . If 
the temperature scale for inflation is assumed to be above 1 TeV, it requires C/A^/^ > 100, 
which for C = 10 implies the self-coupling can be as large as A ~ 0.0001, although it 
also requires a very large number of fields. In any case, it is interesting to examine the 
difficulties that must be overcome in this model to avoid an ultra-flat inflationary potential. 
The fact that such a possibility in remotely realizable is interesting and motivates further 
investigation. 

In summary, the magnitude of A is dependent mutually on the overdamped dynamics 
of warm inflation and the density perturbation requirement. Further development of the 
simple dynamical framework used here may increase A by several orders. 



VII. CONCLUSION 

In this paper, warm inflation solutions have been obtained for the DM model that solve 
the cosmological horizon, flatness and scalar density perturbation problems. Furthermore 
such solution regimes exist for an arbitrarily slow evolution of all macroscopic variables in 
the inflaton field system and the background cosmology. For convenience, the DM models 
considered in this paper had all adjacent mass sites equally spaced, g\Mi+i — Mj| = M. This 
condition is not required to obtain warm inflation solutions. Spacings between adjacent 
sites can widely vary. The inflaton motion remains overdamped for a succession of sites in 
which adjacent spacings are all less than the temperature T. Spacings larger than T are not 
inconceivable, but the inflaton motion is more complicated. For any DM model spacings, 
the ultimate test is the the usefulness of warm inflation solutions that it yields. We have 
examined only the single case of equal spacings. 

There are a few improvements to these calculations that can be made. Our calculations 
have adhered to the high-temperature approximation, which means fields with mass m <T 
are thermally active and those with m > T are thermally dormant. This approximation is 
over restrictive. Generally, fields participate in dissipative and thermalization dynamics once 
m < lOT. An example that treated dissipative heat bath fields with m < 2.5T is in |2^. It is 
worthwhile to extend the present calculation beyond the high-temperature approximation. 
Another improvement to this calculation is to compute the ladder resummed dissipative 



function similar to the shear viscosity case computed by Jeon in [^. This point was noted 
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m 



161 and recently has been verified by Jeon |^ 



The model studied in this paper has been developed into a string theory warm inflation 
p9|,E0[. The DM model has an essential feature for this interpretation. Its 



scenario m 



hierarchy of mass sites are reminiscent of the tower of mass levels of a string. 

The flrst step towards this interpretation was to obtain the DM model from a SUSY 
superpotential ||3^. This has a relevance independent of the string interpretation. It es- 
tablishes that the DM model is natural in the sense of nonrenormalization theorems, which 
means once the parameters are chosen, they stay flxed until SUSY is broken. 

The string picture developed in |3D|,?D| interprets the DM model as an effective SUSY 
model in which the inflaton is a string zero mode and it interacts with higher string mass 
levels, which are the dissipative heat bath flelds. Since the multiplicity of degenerate string 
states increases exponentially with excitation level, strings can provide an adequate supply 
of dissipative heat bath flelds. Interestingly, the dispartity in scales that generally arises 
in DM model warm inflation realizations, Vin^ ^ Tbi > M, is readily explained in the 
string interpretation. Vin^ of the DM model corresponds in the string interpretation to 
the mass scale of a string level for the unperturbed string, i.e., when the coupling to (p is 
switched off. The mass splitting scale M corresponds in the string interpretation to a flne 
structure splitting of a initially degenerate string mass level. This flne structure splitting of 
the level arises from symmetry breaking. As temperature drops below the scale of the string 
mass levels, generally degeneracies in the mass levels will be lifted and thereby create flne 
structure splittings. Since for the warm inflation solutions in this paper Tbi ^ immM, the 
conditions are adequate for various perturbations to break the degeneracy of the mass level. 

For example, if a symmetry breaking occurs at scale iminM > vi > Tbi, the states of 
any mass level split by characteristic scale Vi. Thus after symmetry breaking, the states 
of an initially degenerate mass level shift within a width of order Vi. A generic symmetry 
breaking typically will not lift all the degeneracy, so that a mass level with Af degenerate 
states before symmetry breaking splits into T) < M flnely split levels. In this case, the 
flne structure splitting scale is M ~ Vi/V. Examples of symmetry breaking scenarios and 
estimates of A/", P and Vi are given in [^] . 

For the standard type 1, 11, heterotic and bosonic strings, the string scale is Ms ~ 
lO^^GeV. Therefore the string interpretation requires Vin^ ~ 10^''GeV. For the flrst mass 
level, n = 1, of the heterotic string, for example, TV ~ 2 x 10^ states. V depends on the 
speciflc symmetry breakings that occurs. For a GUT motivated example in we found 
V ~ 10^ 

Comparing these estimates to the bounds in Subsect. |V1 A| , Vin^ is at least one order 
of magnitude below Ms- The example in found VijM ^ 10^. Since Tbi < "^i, to satisfy 
the example in requires Tbi/M < 10^. In Subsect. |V1 A| we found Tbi/M ^ 10^ with 
Tbi ~ 10^°~^^GeV so that vi ~ 10^'^~^^GeV. Thus the estimates from the model in this 
paper are somewhat compatible with the string scenario in |^0|. iminM in our model still is 
below the string scale Ms- Vi is within range of the GUT scales. The improvements to this 
calculation discussed above should elevate iminM to Ms and allow greater flexibility in the 
range of Tbi and M. 

This good news comes at the expense of a flnely tuned self-coupling parameter A. The 
nature of this flne tuning problem is different from a similar problem in supercooled in- 
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flation scenarios Pfl. In our model, the dynamical parameter oc \ip1 is quite large 
relative to the scale of the Hubble parameter. From Subsect. VIA we find in our model 
tti^/Hbi ~ 10^""^. In contrast, in supercooled scenarios m,^ ~ H. Furthermore Subsect. 
VI Q demonstrated that A need not be tiny for observationally interesting warm infiation. 
However this required all the scales to be much smaller and it required an extremely large 
number of fields. Both these requirements confiict with the properties of conventional strings. 
These are disappointing features of this model. Nevertheless, there are indications that in 
this warm infiation dynamics, the density perturbation requirements do not mandate a tiny 
A. Dissipative dynamics provides other means for bounding the density perturbation ampli- 
tude. The present model, despite some of the shortcomings, demonstrates the importance 
of the full decohereing dynamics in determining the density perturbations. 

In this paper, particle production has been treated to the minimal extent of its impli- 
cations for strong dissipation and solving the cosmological puzzles. This comes short of 
attempting a detailed modeling of the produced particle spectrum. This direction has some 
interesting possibilities, which will be discussed here. 

An important fact to note is that the heat bath particles produced in this scenario are 
generally very massive, ~ 10^^ — lO^GeV. Recall in this model, as the infiaton amplitude (po 
relaxes to its minimum, the heat bath fields at a given site i become thermally excited when 
their mass Eqs. (^), ^Xik,'^ik °^ \'fo~ ^il f^-Hs below the temperature and otherwise they 
are thermally unexcited. In order to realize this picture, an elementary requirement is that 
the microscopic time scale of thermalization must be faster than the characteristic time scale 
for the mass of the particle to change, as refiected through the thermal-adiabatic condition 
Eq. (p!9[). In general this condition becomes increasingly less valid for a given heat bath 
field, as its mass becomes increasingly larger than the temperature, because its decay width 
r^^^ decreases exponentially oc exp(— .^/T). Thus, above some m^^^ > T the thermal- 
adiabatic condition no longer will be satisfied. At this point, the existing abundance of the 
given heat bath particles then will freeze-out. After freeze-out, the mass of these heat bath 
particles continues to grow, due to their ipo dependence, but their total number is fixed. Since 
prior to freeze-out the thermal distribution of the heat bath fields also falls exponentially 
oc exp(— m^^^^/T), at time of freeze-out, their abundance could be large if m^^^ ~ T or 
small if m^^jij ~ lOT. As such, the magnitude of the freeze-out abundance depends on how 
quickly the thermal adiabatic condition becomes invalid, and our solutions cover a wide 
range of possibilities. The main point to observe here is that this dynamics generically 
produces very massive particles, in the range ~ 10^^ — lO^GeV once ipo equilibriates. Such 
heavy particle production at the early stages of warm infiation is not important in the post- 
infiation universe, since they will completely dilute during the infiationary period. However, 
this heavy particle productions from the last few mass sites, just before warm infiation ends, 
will leave significant abundances in the post-infiation universe. 

Thus, without further modifications to this model, particle production yields weakly 
interacting massive particles (WIMPs), which are a category of particles that play an im- 
portant role in present-day notions about dark matter [Q. In fact, exceptionally high mass 
particles ~ 10^^ — lO^^GeV, which corresponds to the upper mass range for the heat bath 



fields, sometimes are further distinguished in the literature as "Wimpzillas" |^| 



In general, thermal production of very massive WIMPs is understood to lead to over- 
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abundance problems, which therefore restrict the upper bound on their mass to be ~ lO^GeV 
| 56| . This imphes thermal Wimpzilla production is prohibited for any generic situation, in- 
cluding the model in this paper. Thus most of the mass range for the heat bath fields in our 
model is in danger of an overabundance problem. In fact, as mentioned above, for the model 
in this paper the problem is worse since upon production, the mass of the massive heat bath 
particles can grow by as much as several orders of magnitude due to the y^o-dependence in 
Eqs. (§,(§). 

There are simple remedies to this overabundance problem, which provide a rich number 
of phenomenological possibilities in which the abundance of WIMPs and Wimpzillas can 
be controlled. The simplest solution is to couple the heat bath fields to secondary fields, 
such as the example of a fields in Subsect. VC. This provides decay channels for the 
heat bath fields into less massive or even light fields, thus providing a means to control the 
relative fractions of WIMPs, Wimpzillas, light particles etc... Furthermore, as discussed in 
Subsect. VC, this mechanism can be implemented with minimal effect on the dissipative 
dynamics. Another possibility, which can be applied separately or in conjunction with this, 
is to have a second and short stage of inflation after the initial warm inflation stage, which 
dilutes the WIMP and Wimpzilla abundances. Finally, the lower end of the mass range of 
the heat bath fields found from our calculations falls below the upper acceptability limit for 
thermal WIMP production. The present model has a narrow mass window that falls within 
this acceptability limit. However, with the improvements discussed earlier in this section, 
to the dissipative dynamics, this mass range could be substantially enlarged. 

The calculations in this paper were guided by observation and consistency with quantum 
field theory. "Nice" particle physics was not an a priori requirement. On the one hand, this 
is a constructive approach that attempts to find a toy model to study a very complicated 
dynamics. On the other hand, this is a predictive approach. Quantum field theory is 
postulated to hold at the scales of inflationary dynamics. Combining this theoretical tool 
with observation, a mutually consistent model has been deduced. 

We can conclude that this model is a good constructive tool for studying warm inflation 
dynamics. No complete conclusion is possible as yet about the predictive content of this 
model. The interesting connection between this model and strings found in [p9| , ^ implies a 
revised meaning is required of "nice" particle physics. At the inflationary scale, nice particle 
physics may not be synonymous with simple particle physics. This model is evidence that 
inflationary dynamics can be a multi- or even ultramulti- field problem. 
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FIGURE CAPTION 



Figure 1: Warm inflation scales as a function of the adiabatic parameter a for solid 
lines: iV = 5, A^^ = 65, /? = 0.5 and C = 1/8; dashed lines: A^^^ = 1/8, A^^ = 65, /3 = 0.5 
and C = 1/8. For the scalar density perturbations, 5p/ p = 10~^ (< 10^^) to the left (right) 
of the vertical solid and dashed lines for the two corresponding cases. For the entire a range, 
the spectrum of density perturbations is flat up to logarithmic corrections. 
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